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ABSTRACT 


Ice  covers  cause  a number  of  changes  in  alluvial  channel  flows 
by  approximately  doubling  the  wetted  perimeter  and  thereby  producing  a 
redistribution  of  the  boundary  and  internal  shear  stresses.  A series 
of  flume  experiments  was  conducted  to  investigate  the  effects  of  simu- 
lated ice  covers  on  the  following  characteristics  of  alluvial  channel 
flows:  (1)  depth  versus  discharge  and  friction  factor  relationships; 

(2)  velocity  distribution;  (3)  suspended  sediment  concentration  profiles 
and  discharge;  (4)  total  bed-material  discharge;  and  (5)  characteristics 
of  bed  forms  and  their  relationship  to  bed  friction  factors. 

In  comparison  to  free-surface  flows  with  the  same  unit  dis- 
charge and  energy  slope,  flows  with  simulated  ice  covers  were  found  to 
have  substantially  larger  depths  and  lower  average  velocities.  Due 
mainly  to  the  lower  velocities,  sediment  discharges  were  found  to  be 
sharply  reduced.  The  flow  in  an  ice-covered  channel  is  divided  by  a 
plane  of  zero  shear  stress  (locus  of  maximum  velocities)  into  a lower 
and  an  upper  layer,  each  with  average  velocities  approximately  equal  to 
the  overall  average  velocity.  To  the  extent  that  the  shear-stress  and 
velocity  distributions  in  the  lower  layer  of  an  ice-covered  flow  are  the 
same  as  in  a free  surface  flow  with  the  same  mean  velocity, and  depth 
equal  to  the  thickness  of  the  lower  layer,  relationships  developed  for 
flows  in  alluvial  channels  with  a free  surface  can  be  used  for  predic- 
ting certain  features  of  flows  in  ice-covered  alluvial  channels.  This 
applies  in  particular  to  bed  friction  factors,  sediment  discharge,  and 
apparently  also  to  bed-form  characteristics. 


I.  INTRODUCTION 


A.  Background.  At  the  higher  latitudes  in  Earth's  northern 
and  southern  hemispheres,  the  surfaces  of  most  natural  bodies  of  water 
undergo  an  annual  period  of  freezing.  In  alluvial  channels,  the  pres- 
ence of  an  ice  cover  radically  alters  many  features  of  the  water-  and 
sediment-transport  processes.  The  dynamics  of  this  type  of  system  are 
none  too  well  understood  even  for  the  simpler  case  of  flow  in  an  ice- 
free  alluvial  channel,  mainly  because  it  is  the  interaction  between 
the  flow  and  the  channel  bed  that  determines  the  configuration  of  the 
bed  and  consequently  its  hydraulic  roughness  and  the  depth  of  flow. 

Thus  in  contrast  to  a rigid-boundary  channel,  wherein  only  the  wa:er 
surface  is  a free  surface,  an  alluvial  channel  has  two  free  surfaces.  In 
an  ice-covered  alluvial  channel,  this  complexity  is  compounded  further 
by  the  variability  of  the  hydraulic  roughness  of  the  underside  of  the  ice 
cover.  The  under  surface  is  typically  very  smooth  during  formation  and 
growth  of  the  ice  cover  when  temperatures  are  cold,  but  becomes  rippled 
following  the  onset  of  higher  temperatures  and  melting  from  below,  and 
becomes  very  rough  after  candling  and  break-up  of  the  cover  occur  (Ashton 
and  Kennedy  1972) . Thus  not  only  the  presence  of  an  ice  cover,  but  also 
any  change  in  the  roughness  of  its  underside  is  expected  to  cause  changes 
in  the  bed  configuration  of  the  channel  and  in  the  flow  and  sediment 
transport  characteristics  of  the  stream. 

Some  of  the  qualitative  effects  of  a stationary,  floating  ice 
cover  on  the  flow  and  sediment  transport  regimes  can  be  inferred  from 
elementary  physical  reasoning.  The  addition  of  a solid  upper  boundary 
surface  approximately  doubles  the  wetted  perimeter  of  the  channel.  Bar- 
ring a compensating  reduction  in  the  hydraulic  roughness  of  the  bed, 
the  total  boundary  shear  and  consequently  the  depth  of  flow  will  be 
increased,  and  the  average  velocity  will  be  reduced.  Because  the  total 
boundary  shear  must  be  divided  between  the  upper  and  lower  surfaces,  a 
reduction  in  shear  stress  at  the  bed  is  to  be  expected  unless  the  depth 
is  roughly  doubled.  The  reduction  in  velocity  and  bed  shear  stress  is 
in  turn  expected  to  reduce  both  the  bed  and  the  suspended  sediment  loads. 


Changes  in  the  internal  distributions  of  turbulent  shear  stress  and 
velocity  are  apt  to  be  even  more  drastic.  The  redistribution  of  turbulent 
shear  stress  will  affect  the  turbulent  mass  diffusivity  and  consequently 
the  vertical  distribution  of  suspended  and  probably  also  the  ratio  of 
suspended  load  to  bed  load. 

Quantitative  prediction  of  the  ice-cover  effects  would  require 
knowledge  of  the  magnitudes  and  distributions  of  the  new  boundary  and 
internal  turbulent  shear  stresses,  which  in  turn  depend  on  the  relative 
degrees  of  hydraulic  roughness  of  the  underside  of  the  ice  cover  and 
the  channel  bed.  The  principal  difficulty  lies  in  the  freedom  of  the 
bed  configuration  to  adjust  to  the  new  shear  stress  distribution,  which 
in  turn  is  modified  by  changes  in  the  bed  configuration.  in  the  parlance 
of  system  dynamics,  the  stream  contains  an  internal  feedback  loop. 

Progress  toward  development  of  a quantitative  model  for  predicting  the 
effects  of  an  ice  cover  on  alluvial  channel  flow  processes  requires  a 
better  understanding  of  how  this  feedback  loop  operates. 

The  river-morphology  literature  contains  very  little 
information  on  ice-cover  effects  that  could  serve  as  a guide  in  develop- 
ing a quantitative  predictive  model.  In  addition,  the  accuracy  of  the 
limited  amount  of  sediment-transport  data  for  ice-covered  streams  that 
is  available,  is  open  to  question  because  of:  (1)  the  reduced  relia- 
bility of  sampling  equipment  and  methods  when  operating  under  freezing 
conditions;  and  (2)  insufficient  knowledge  concerning  the  interaction 
between  the  ice  cover  and  the  flow  to  interpret  fully  the  measurements 
and  design  an  efficient  and  reliable  sampling  program.  Heretofore,  these 
problems  have  been  investigated  little,  both  because  of  the  many  dif- 
ficulties entailed  in  obtaining  sediment  transport  and  related  measure- 
ments under  freezing  conditions,  and  a general  consensus  that  flow  dis- 
charges and  sediment  transport  rates  are  much  reduced  when  ice  covers 
exist . 

The  motivation  for  the  research  described  in  this  report  is 
an  outgrowth  of  the  U.S.  Geological  Survey's  charge  to  develop  the 
National  Stream  Quality  Accounting  Network  (NASQUAN).  To  meet  the 
NASQUAN  objective  of  describing  the  water  quality  of  the  Nution's  streams 
and  rivers  on  a systematic  and  continuing  basis,  winter  data  acquisition 


and  interpretation  activities  in  parts  of  the  country  with  a freezing 
climate  need  to  be  upgraded  and  expanded. 
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B.  Purpose  and  Scope  of  Investigation.  The  main  purpose 
of  this  investigation  was  to  obtain  basic  information  on  how  ice 
covers  affect  the  sediment-transport  characteristics  of  alluvial 
streams.  As  outlined  in  section  I-A,  these  characteristics  include  the 
mean  velocity  and  depth  of  flow,  the  velocity  distribution,  the  hydraulic 
roughness  of  the  bed  and  ice  cover,  the  magnitude  and  vertical  distribu- 
tion of  the  suspended  sediment  load,  and  either  the  bed  or  the  total 
sediment  load.  The  various  ice-cover  effects  are  investigated  and  analyzed 
with  a view  toward  formulating  recommendations  that  would  help  to  improve 
field  sampling  and  data  collection  programs  for  ice-covered  alluvial 
channels. 

The  project  activities  have  included:  (1)  a one-day  USGS-I1HR 
Workshop  on  Measurement  of  Suspended  Sediment  Transport  in  Ice-Covered 
Alluvial  Streams  held  at  The  University  of  Iowa  on  November  25,  1975; 

(2)  a field-trip  to  the  Upper  Iowa  River  near  Dorchester,  Iowa  on  January 
20,  1976,  to  observe  a USGS  flow  and  sediment  discharge  measurement  in 
an  ice-covered  stream;  and  (3)  a series  of  experiments  in  a laboratory 
flume  wherein  detailed  comparisons  were  made  between  the  characteristics 
of  alluvial-channel  flow  processes  occurring  under  simulated  ice  covers, 
and  those  occurring  in  equivalent  free  surface  flows. 

Aside  from  brief  reviews  in  the  Appendices,  this  report  is 
restricted  to  the  laboratory-flume  experimental  investigation.  The 
experimental  investigation  addresses  the  effects  of  smooth  and  rough 
ice  covers  on  the  following  characteristics  of  the  flows:  (1)  depth 
versus  discharge  and  friction-factor  relationships;  (2)  velocity  dis- 
tribution; (3)  suspended  sediment  concentration  distributions  and 
discharge;  (4)  total  bed-material  discharge;  and  (5)  characteristics 
of  bed  forms  and  their  relationship  to  bed  friction  factors. 


II.  BASIC  HYPOTHESES 


Throughout  this  investigation,  flows  occurring  under  an 
ice  cover  are  compared  insofar  as  possible  with  equivalent  free-sur- 
face  flows.  The  concept  of  equivalent  ice-covered  and  free-surface 
flows  as  used  herein  embodies  pairs  of  two-dimensional  uniform  flows, 
one  with  a free  surface  and  the  other  with  an  ice  cover,  for  which  the 
unit  water  discharge,  q,  the  slope,  S,  and  the  fluid  and  sediment  prop- 
erties are  all  the  same,  but  the  mean  depth,  d,  hydraulic  radius,  r, 
mean  velocity,  U,  friction  factor,  f (or  Manning's  n) , and  the  unit 
total  sediment  discharge,  qgt»  are  different.  The  ice  cover  is  assumed 
to  be  stationary  and  floating  so  that  the  pressure  at  the  ice-water 
interface  is  equal  to  the  atmospheric  pressure  plus  the  weight  of  an 
ice  column  of  unit  area  and  height  equal  to  the  thickness  of  the  ice 

cover.  In  nature,  this  corresponds  to  a comparison  between  the  state  of 

affairs  which  would  exist  prior  to  and  for  some  time  following  the  forma- 
tion of  an  ice  cover  in  a long  reach  of  channel  with  steep  banks  (i.e., 
no  appreciable  increase  in  channel  width  due  to  the  increased  depth  caused 

by  the  ice  cover).  For  both  flows  q is  assumed  to  be  a dependent 

variable,  and  longer  term  effects  due  to  changes  in  the  imposed 
sediment  load,  or  changes  in  slope  or  flow  distribution  due  to  shifts 
in  meander  patterns,  for  example,  are  not  considered.  Other  phenomena 
relating  to  sediment  transport  which  are  not  considered  include  ice  jams 
and  the  uptake  and  transport  of  bottom  sediments  by  clusters  of  frazil 
ice  and  dislodged  chunks  of  anchor  ice. 


A.  Depth  of  Flow.  The  Chezy  resistance  equation  for  a uni- 
form flow  in  a wide  channel  may  be  expressed  as 


q 


= C d/rs 


(1) 


in  which  the  Chezy  coefficient,  C, 
C - ! r1/6 


is 


n 

according  to  the  Manning  formula  (in  metric  units),  where  n is  the 
Manning  friction  factor,  and 


(2) 
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C - /8g/r  (3) 

according  to  the  Darcy-Weishach  formula,  where  f is  the  Darcy-Welsbach 
friction  factor,  and  g is  the  acceleration  of  gravity.  Using  the  sub- 
script I to  denote  an  ice-covered  flow,  and  the  subscript  F to  denote 
an  equivalent  free-surface  flow,  the  hydraulic  radii  for  the  respective 
flows  are  r^  - d/2  and  r^  - d.  Since  * q^  and  Sj  * S^,  for  equivalent 
flows,  the  foregoing  relationships  can  be  combined  to  give  the  ratio 
of  the  depth  of  the  ice-covered  flow  to  the  depth  of  the  equivalent  free- 
surface  flow,  which  is 


. 20'1  (V-6  • 1.32  A)0'6 

dF  nF  nF 


according  to  the  Manning  formula,  and 

f 


. o 1 / 3 /lxl/3 

d 'f  ' 

F F 


1.26  (-p)1/3 


(4) 


(5) 


according  to  the  Darcy-Weisbach  formula.  In  (4)  and  (5)  n^  and  f^ 
are  composite  friction  factors,  reflecting  the  roughness  of  both  the 
bed  and  the  under  surface  of  the  ice  cover.  Equations  4 and  5 indicate 
that  unless  the  composite  friction  factors  for  an  ice-covered  flow 
are  substantially  smaller  than  the  bed  friction  factors  for  an  equiva- 
lent free-surface  flow,  the  depth  of  flow  will  be  significantly  in- 
creased by  the  formation  of  an  ice  cover.  Assuming  the  same  bed  rough- 
ness for  ice-covered  and  free-surface  conditions,  Carey's  (1966)  measure- 
ments in  the  St.  Croix  River,  Wisconsin,  indicated  values  of  n /n  ranging 
from  0.73  to  1.06  which  would  result  in  di/dp  ranging  from  1.09  to 
1.37. 

Discussion  pertaining  to  how  the  composite  friction  factors 
n^  and  f^  are  related  to  the  friction  factors  for  the  bed  and  the  under- 
surface of  the  ice  cover  is  deferred,  pending  discussion  of  the  effect 
of  an  ice  cover  on  the  vertical  distribution  of  velocity. 


' — ! 
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B.  Velocity  Distribution.  The  Karman-Prandtl  velocity  de- 
fect law  for  turbulent  flow  written  in  the  form 


U*i  yi 

Ui  Umax  tc^  Y^ 


(6) 


together  with  the  equation  for  the  linear  vertical  distribution  of 
shear  stress 


(7) 


provide  an  appropriate  model  for  the  internal  structure  of  a uniform 
two-dimensional  flow  in  an  ice-covered  channel.  Most  of  the  quantities 
in  (6)  and  (7)  are  represented  in  figure  1,  which  illustrates  the  verti- 
cal distribution  of  velocity,  u^,  and  shear  stress  x^.  The  subscript 
i ■ 1,  2 denote  the  lower  or  upper  layer,  respectively.  Previously 
undefined  quantities  in  (6)  and  (7)  are:  the  maximum  velocity,  umax» 
which  occurs  at  the  Interface  between  two  layers,  where  x^  “ 0;  the 
distance,  y^,  from  the  lower  or  upper  boundary;  the  thickness,  Y^,  of 
the  lower  or  upper  layer;  the  von  Karman  turbulence  coefficient,  k^; 
the  specific  weight  of  water,  y;  and  the  shear  velocity 


u*i  - - s#? 

where 


(8) 


T0i  “ yYiS 


(9) 


is  the  shear  stress  at  the  lower  or  upper  boundary,  and  p is  the  fluid 
mass  density.  Also  from  figure  1 and  (9)  it  is  clear  that  the  average 
or  composite  boundary  shear  stress  for  an  ice-covered  flow,  taking  both 
the  lower  and  upper  boundaries  into  account,  is 


T01  + T02 


01 


- ydS/2 


(10) 


Equations  6 and  7 can  both  be  obtained  from  the  x-direction 
Reynolds'  equation  for  steady,  uniform  turbulent  flow  between  two  paral- 
lel plates  that  are  oriented  horizontally  in  the  transverse  direction 


Shear  Stress 


Velocity 


Definition  sketch  for  velocity  and  shear  stress 
distributions  in  an  ice-covered  channel. 


Figure  1 
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and  parallel  to  the  x axis  which  is  tilted  downward  at  the  slope  S, 
„ du,  

- YS 


dy 


— (pv 


dy. 


pu|v|) 


P8, 


(ID 


ri  “'i 

Under  these  conditions  (11)  is  nothing  more  than  the  statement  of  static 
equilibrium  for  a fluid  element  given  by  (7),  where  is  the  sum  of  the 
viscous  and  turbulent  shear  stresses.  Equation  6 can  be  obtained  by 
integrating  (11)  first  from  0 to  y and  then  from  y^  to  Y^,  following 
introduction  of  the  Boussinesq  eddy  diffusivity,  e^,  as  defined  by 


dy. 


- ujv’ 


where 


y4 

•c.u  y (1  - — ) 
1 i 1 


(12) 


(13) 


and  the  boundary  conditions 

du. 


pv 


T0i  " YYis 


V ui 


max 


x^  ■ 0 


Also,  in  the  second  integration  of  (11),  the  viscous  shear  stress  is 
neglected  since  v « except  in  the  immediate  vicinity  of  y^  =0. 
The  mean  velocity  in  the  lower  or  upper  layer  is  given  by 


U 


i 


Uidyi 


(14) 


which  can  be  used  to  evaluate 

Essentially  the  same  model  for  the  internal  flow  structure 
was  used  by  Hancu  (1967),  as  reported  by  Uzuner  (1975),  and  by  Larsen 
(1969)  to  describe  the  flow  in  an  ice-covered  channel,  and  also  by 
Ismail  (1952)  in  a study  of  sediment  suspension  in  a closed  rectangular 
channel.  Clearly  the  model  applies  equally  well  to  a free  surface  flow, 
in  which  case  Y^  * d,  the  total  depth  of  flow,  and  \ ^ * 0. 
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C.  Friction  Factors.  The  relative  magnitudes  of  the  fric- 
tion factors  for  the  bed  and  the  underside  of  the  ice  cover  are  impor- 
tant because  they  influence  the  division  of  shear  stress  between  the 
lower  and  upper  boundaries,  which  in  turn  fixes  the  relative  magnitudes 

of  Y and  Y_  and  the  vertical  position  of  u . The  Manning  (metric 
i i.  max 

system)  and  Darcy-Weisbach  friction  factors  for  the  lower  and  upper 
surfaces  are  related  to  the  bulk  flow  characteristics  in  the  lower  and 
upper  layers  by  the  respective  equations 


,2/3  gl/2 


,5/3  gl/2 


Y Y 

f - 8gS  — ^ - 8gS  —^2 
Ui  «i 

The  corresponding  composite  friction  factors  for  a wide  ice-covered 
channel,  for  which  r - d/2,  are 

(d/2)2/3  S1/2 

n ■ '■*■■■'  ■ • 


fI  - AgS 

U 

In  terms  of  the  friction  factors  for  the  upper  and  lower  surfaces,  the 
composite  friction  factors  given  by  Hancu's  (1967)  equations,  as 
reported  by  Uzuner  (1975),  are 

i 


2,U1.2  d .1/3  2,°2.2  ,d  .1/3  1/2 

"i  ■ (i>  hi(u->  (r>  + n2<ir)  (r> 

LI  2 

and 

fi  ■ 1 |<r)2  fi + <r>2  f2]  <» 

One  of  the  main  questions  to  be  investigated  in  the  present 
study  involves  the  extent  to  which  bed  friction  factors,  characteristics 
of  bed  forms,  and  sediment-transport  relationships  in  an  ice-covered 


i 


channel  are  equivalent  to  those  in  a free-surface  flow  with  depth  and 
mean  velocity  U^.  This  is  a reasonable  hypothesis  insofar  as  y^  * 
is  a plane  of  zero  shear  stress  like  the  water  surface  in  a free-surface 
flow.  The  relationships  of  Yj  and  to  the  velocity  and  shear-stress 
distributions  and  bed  friction  factor  in  the  lower  layer  of  an  ice- 
covered  flow  then  would  be  exactly  the  same  as  those  of  d and  U to  the 
corresponding  features  of  a free  surface  flow.  It  is  therefore  important 
to  be  able  to  estimate  Y^/d,  the  fractional  depth  where  maximum  velocity 
occurs,  in  an  ice-covered  flow.  If  the  friction  factors  for  both  sur- 
faces can  be  estimated,  and  q,  S,  and  d are  known,  this  can  be  done 
by  solving  the  equations 


YL/d  + Y2/d  - 1 


(21) 


and 


qx  + q2  * q 


(22) 


simultaneously  by  trial  and  error  following  the  substitution  of  either 
(15)  or  (16)  into  (22)  to  obtain  either 


(Yx/d)5/3  (Y2/d)5/3 

nl  n2  d5^3  S1^2 


(23) 


or 


(Y  /d)3/2  (Y  /d)3/2 

« 1/2  + , 1/2 


'8gS  d 


3/2 


(24) 


If,  as  a first  approximation,  it  is  assumed  that  - U2  * U,  the  same 
procedure  results  in  the  simpler  approximate  equations 

Y, 


or 


r:u  + (^>3/2i_1 


(25) 


fl  +f2 


(26) 
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D.  Suspended  Sediment.  For  a steady,  uniform  flow  in  an 
ice-covered  channel,  as  for  a free-surface  flow,  the  vertical  distri- 
bution of  suspended  sediment  is  governed  by  the  Schmidt-0 ' Brien  equa- 
t ion 

w C + c ^ - 0 (27) 

p s dy 


which  expresses  a balance  between  downward  transport  due  to  gravita- 
tional settling  and  net  upward  transport  due  to  turbulent  diffusion. 

In  (27),  C is  the  concentration  of  suspended  sediment  (usually  expressed 
in  mg/fc),  Wp  is  the  terminal  settling  velocity  of  the  sediment  particles, 

and  e is  the  vertical  turbulent  mass  transfer  coefficient  for  the  sus- 
s 

pended  sediment.  It  is  assumed  that  eg  is  related  to  the  eddy  diffusivity 

e “ e,  in  (13)  by  the  relation 
mo  i 


e 

s 


6e_ 


(28) 


where  8 is  the  reciprocal  of  the  turbulent  Schmidt  number.  For  fine 
sediment  particles  (w^/ku*  < 0.3),  it  is  usually  accepted  that  8 1 1. 

For  larger  particles  there  is  conflicting  evidence  (e.g.,  Jobson  and 
Sayre, 1970)  as  to  whether  8 is  smaller  or  larger  than  one,  or  oven 
whether  it  is  constant  over  the  depth  of  flow. 

The  model  for  the  vertical  distribution  of  c_  adopted  herein 
for  a flow  with  an  ice  cover  is  an  extension  of  the  one  developed  by  Ismail 
(1952)  in  an  experimental  investigation  of  suspended  sediment  transport 
in  a closed  rectangular  conduit.  A definition  sketch  is  shown  in  figure 
2,  where  n » y/d  is  the  normalized  elevation  above  the  bed.  Near  the 
lower  and  upper  boundaries  cg  follows  the  parabolic  distribution  usu- 
ally associated  with  a logaritlunic  velocity  distribution.  However,  in 
the  mid-depth  region,  cg  is  assumed  to  be  distributed  linearly  between 
the  vertices  of  the  lower  and  upper  parabolas.  The  model  depicted  in 
figure  2 differs  from  Ismail's  in  two  respects.  First,  the  one  shown 
in  figure  2 allows  for  different  upper-  and  lower-boundary  roughnesses 
so  that  symmetry  about  n ” 0.50  is  not  required.  Second,  the  constant  value 
of  e^  adopted  by  Ismail  for  the  mid-depth  region  is  twelve  percent  less 
than  maximum  value  at  the  vertices. 


coefficient  for  suspended  sediment  in  an  ice-covered  channel. 


(35) 


0(1  - Vd>  2 


where  n.  - 1 - Vj/d.  z2  - 

L Z.  *2 


The  depth-averaged  value  of  e , 

s 

over  the  three  layers,  is 


obtained  by  integrating 


e 

s 


e 

s 


1 

/ e dn 

o s 


Vlu* 


Y1  1 
Y (— =—  + — ) 

1 Vl2d  16; 


+ 62"2“.,  Y2(T2d  + y <*> 


For  a free-surface  flow,  (30)  and  (31),  with  Y^  » <i  apply  over 
the  whole  range  0 < n < 1 . The  depth-averaged  value  of  eg  for  this 

case  is 

e * 6 tcu^d/6  (37) 

Equations  31,  33,  and  35  all  express  the  vertical  distribution 

of  C relative  to  the  concent  rat  Ion  C . often  determined  bv  measurement. 

a * 

at  some  reference  level  n » Lack  of  a reliable  method  tor  predicting 

C corresponding  to  a given  set  of  t . >w  atu:  ,.u  i.  ter  l sties  is 

a 

one  of  the  weaknesses  in  current  sediment  t m-s  r,.  ties.  Einstein 

(1950)  attempted  to  relate  C at  .»  -list  u.  . ! . luesters 

a 

above  the  bed  to  the  theoretical  tee  . . ...*.  it.  .it  the 

outer  edge  of  the  viscous  sublayer.  Wi  n.ivt  pro- 
posed relating  C to  the  theoreli.al  ...  ...  •.  ..-.int  tlux, 

a 

the  concentration  in  the  bed,  and  th«  probat  ...t.  : t.  i;  ..  . , al  bed  ele- 
vation equals  or  exceeds  the  elevation  above  to  w:  i ev»  . . v , at 

' m 

which  the  maximum  flux  occurs. 

Using  Ismail's  (1952)  model  for  , Shell  and  Harden  (1978) 
predicted  profiles  with  lower  suspended  sediment  concentrations  tor  <ce- 
covered  than  for  free-surface  flows,  especially  in  the  upper  region. 

They  attributed  this  to  a reduction  in  the  vertical  mixing  capacity  due 
to  the  presence  of  the  ice  cover. 


E.  Total  Bed-Material  Discharge.  When  computed  for  a parti- 
cular channel,  unit  discharges  of  total  bed-material  load,  q t>  and 
bed  load,  q^*  predicted  by  most  of  the  currently-used  formulas,  can  be 


..  ’ 


approximated  by  a generalized  relationship  of  the  type 


C „ } “ <U  - V”  <T0  - V"  < 

sb 

The  exponents  m and  n have  different  values  (usually  between  0 and 

about  3,  but  with  m+n  * 2 to  3)  in  the  different  formulas,  and  U and 

c 

tc>  if  included,  represent  the  critical  mean  velocity  and  bed  shear 
stress  at  which  particles  begin  to  move. 

For  the  equivalent  ice-covered  and  free-surface  flows  as 
defined  at  the  beginning  of  this  chapter  the  mean  velocity  in  the  ice- 
covered  flow  will  be  reduced  according  to  the  ratio 


and  the  bed  shear  stress  according  to  the  ratio 


Y1  d] 
= (— ) — 
(d  ’l  dT 


With  typical  values  for  both  ratios  ranging  between  about  0.7  and  0.9, 
it  is  clear  from  (38)  that  bed-material  and  bed-load  transport  rates 
can  be  greatly  reduced  by  the  presence  of  an  ice  cover. 

Fundamental  work-energy  considerations  lead  to  the  same 
conclusion.  The  equation  of  motion  for  steady  uniform  flow,  integrated 
over  the  cross-sectional  area,  can  be  reduced  to  the  statement 

yQS  = U 7^  P (4 

where  the  left-hand  side  represents  the  mean  rate  of  energy  dissipation 
per  unit  length  of  channel,  and  the  right-hand  side  represents  the  mean 
rate  of  work  to  overcome  resistance  at  the  boundary,  denoted  by  the 
wetted  perimeter,  P.  Considering  equivalent  ice-covered  and  free-sur- 
face flows,  the  left-hand  side  is  the  same  for  both  cases.  However 
the  right-hand  side  for  the  free-surface  flow  is 


RHSF  = UF  T0F  PF 


where  the  wetted  perimeter,  P^.,  consists  of  the  bed  (including  submerged 
banks),  whereas  for  the  ice-covered  flow  it  is 


RHS, 


U1  T01  P1  + U1  T02  P2 


(43) 


where  the  work  is  divided  between  the  wetted  perimeter  of  the  bed,  P 1 p , 
and  the  wetted  perimeter  of  the  ice  cover,  P^.  Since  the  two  right-hand 
sides  have  the  same  value,  it  is  clear  that  the  energy  dissipated  by 
the  bed  forces  must  be  significantly  less  in  the  ice-covered  flow.  The 
fraction  of  the  total  work  expended  in  overcoming  bed  resistance  is  equal 
to  X01  P1^T01  P1  + t02  P2^’  which  reduces  to  Y^/d  for  a wide  rectangular 


channel  where  P 


1 


P2*  T01  = YY1S’  and  t02  = yY2S‘ 


It  is  reasonable  to  assume  chat  q . and  the  entrainmenc  of  bed 

sb 

material  into  suspension  are  associated  only  with  that  fraction  of  the 
work  that  is  exerted  in  overcoming  bed  resistance.  It  follows  that  both 
aud  the  suspended  sediment  unit  discharge 


qss  = /. 


0+ 


uC  dy 


(44) 


must  be  significantly  reduced  by  the  presence  of  an  ice  cover.  With 
respect  to  (44),  one  would  not  expect  the  increase  in  d,  due  to  the  ice 
cover,  to  be  sufficient  to  offset  the  reduction  in  C,  due  to  the  decreased 
capacity  of  the  flow  to  entrain  sediment  from  the  bed  and  maintain  it 
in  suspension. 

The  foregoing  conclusions  concerning  the  reduction  in  bed- 
material  transport  due  to  the  presence  of  an  ice  cover  are  consistent 
with  observations  on  a number  of  rivers  in  Canada  reported  by  Tywoniuk 
and  Fowler  (1972).  They  found  both  flow  and  suspended  sediment  trans- 
port to  be  usually  much  lower  and  less  variable  during  the  ice-cover 
periods,  with  roughly  10  percent  of  the  total  annual  suspended  sediment 
discharge  transported  by  the  20  percent  of  the  total  annual  flow  occurring 
during  the  approximately  5-month  ice-cover  period.  Observations  which 
contradicted  this  trend  were  obtained  mainly  in  break-up  periods, 
during  which  concentrations  and  discharges  may  be  relatively  high. 


F.  Bod  Forms.  In  a given  sand  bed  channel,  characterized  by 
a given  slope,  width,  and  distribution  of  bed-material  particle  sizes, 
the  bed  forms  progress  through  an  orderly  spectrum  of  configurations, 
from  ripples  to  dunes  to  flat  bed  to  antidunes  to  chutes  and  pools  as 
the  water  and  sediment  discharges  are  increased.  The  heights  and  lengths 
of  these  bed  forms  are  responsive  to  changes  in  the  flow  conditions,  and 
in  turn  the  resulting  bed  configuration  largely  determines  the  hydraulic 
roughness,  thereby  significantly  influencing  the  depth  and  velocity  of 
flow.  The  way  in  which  the  presence  of  an  ice  cover  affects  the  bed 
forms  is  not  clear  from  presently  available  bed-form  predictors,  all  of 
which  have  been  developed  for  free-surface  flows. 

Vanoni  (1974)  developed  a set  of  charts,  each  for  a different 
range  of  bed-material  sizes,  which  consist  of  graphs  of  the  Froude  number, 
IF-  U/tgd,  versus  the  ratio  of  depth  to  median  grain  diameter,  d/D^,  on 
which  the  data  group  into  zones  corresponding  to  different  types  of 
bed  forms.  In  general  the  bed  forms  progress  toward  the  upper  end  of 
the  bed-form  spectrum  from  ripples  to  dunes  to  transition  to  flat  bed 
to  antidunes  to  chutes  and  pools  as  IF  or  d/D^  increases.  In  using  the 
charts  for  ice-covered  channels,  it  is  not  entirely  clear  whether  d or 
Y(  is  the  more  appropriate  length  scale  for  the  flow  in  IF  and  d/Il,^, 
or  even  whether  the  same  length  scale  should  be  used  in  both  parameters. 

In  particular,  the  physical  significance  of  IF’  is  not  clear  for  an  Ice- 
covered  flow.  If  IF  is  interpreted  as  ► 8S / f where  f - f , then  it  ap- 
pears reasonable  to  use  Y ^ . If  Y^  is  used  in  both  Fund  d/l),^,  forma- 
tion of  an  ice  cover  will  cause  both  F and  Y^/D^  to  decrease  below 
their  corresponding  values  for  an  equivalent  free-surface  flow.  Indi- 
cating a regression  toward  the  lower  end  of  the  bed-form  spectrum.  If 
d^  is  used  in  both  parameters,  Fwlll  decrease  but  d/I),  ^ will  increase, 
although  probably  not  by  enough  to  offset  the  decrease  in  F,  so  an  indi- 
cation of  regression  toward  the  lower  end  of  the  bed-form  spectrum  is 
still  likely.  However,  if  Y^  is  used  in  F and  d^  in  d/1)^^,  progression 
toward  the  upper  end  of  the  spectrum  would  bo  indicated  in  some  cases. 

Simons  and  Richardson  (1966)  presented  a graphical  relationship 
of  stream  power,  i^U,  versus  median  particle  fall  diameter  wherein,  for 
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any  given  fall  diameter,  the  progression  toward  the  upper  end  of  the  bed- 
form  spectrum  is  always  in  the  direction  of  increasing  stream  power.  In 
light  of  the  discussion  on  work-energy  considerations  in  the  last  sec- 
tion it  seems  clear  that  only  that  part  of  the  stream  power  that  is 
associated  with  the  bed  (i.e.,  gY^SU)  should  be  considered  when  using 
the  relationship  for  ice-covered  flows.  Because  the  stream  power  asso- 
ciated with  the  bed  is  always  less  for  an  ice-covered  flow  than  for  an 
equivalent  free-surface  flow,  regression  toward  the  lower  end  of  the 
bed-form  spectrum  accompanying  formation  of  an  ice  cover  would  invariably 
be  predicted  by  the  Slmons-Rlchardson  relationship. 

Willis  and  Kennedy  (1977)  presented  graphical  relationships 

wherein  the  ratios  of  indicators,  obtained  by  spectral  analysis,  of 

mean  length,  L,  and  mean  height,  H,  of  bed  forms  to  depth  are  shown  to 
2 

increase  with  11’  in  the  sub-transition  range  of  the  bed-form  spectrum. 

In  applying  these  relationships  to  ice-covered  flows,  it  is  again  un- 
certain whether  Y^  or  d^  is  the  more  appropriate  length  scale  to  repre- 
sent the  flow.  Calculations  using  typical  values  of  d /d  and  Y /d  show 

Jr  lr 

that  the  relationships  would  predict  a reduction  in  both  length  and 
height  of  bed  forms  accompanying  format 'on  of  an  ice  cover  for  all  com- 
binations of  Yj  and  d^  in  IF"  and  (1,  or  H)/d,  except  when  d^  is  used  in 
the  former  and  Y^  in  the  latter.  However  there  is  considerable  variation 
in  the  results  predicted  using  the  other  three  combinations. 

Grigg  (1970)  reported  data  for  free-surface  flows  with  ripples 
and  dunes.  Indicating  that  both  L and  o (standard  deviation  of  hed 

y 

elevation,  which  is  roughly  proportional  to  H)  increase  with  stream  power. 
For  a given  grain  size,  the  data  also  indicate  a fairly  consistent 
tendency  for  L/o^  to  increase  with  stream  power.  Considering  only  that 
portion  of  the  stream  power  that  is  expended  in  overcoming  bed  resistance, 
one  would  expect  from  these  trends  that  formation  of  an  ice  cover 
would  invariably  cause  a reduction  in  the  size  of  the  bed  forms,  but 
possibly  cause  the  hydraulic  roughness  of  the  bed  to  increase  due  to 
a reduction  in  L/H. 

Alam  and  Kennedy  (1969),  in  an  Investigation  of  friction 
factors  for  flow  in  sand  bed  channels,  used  dimensional  analysis  together 
with  flume  and  river  data  obtained  by  several  investigators  to  develop 
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a chart  wherein  the  bed-form  friction  factor,  f",  is  represented  as 

a function  of  ll//eb_,,  and  r./Dr,,,  where  r is  the  hydraulic  radius 
" 50  b 50  b 

associated  with  the  bed.  Replacing  f"  by  f"  and  r^  by  Vj  for  flows 
with  ice  covers,  the  chart  indicates  that  in  most  cases  i j'  would  be 
larger  than  f"  for  an  equivalent  free-surface  flow.  The  Increase  in 
friction  factor  could  be  due  to  a reduction  in  L/H  even  through  there 
is  a decrease  in  both  L and  H. 

C. ■ Summary.  For  easy  reference,  the  main  items  introduced  in 
this  chapter  that  pertain  to  the  effects  of  an  ice  cover  on  alluvial 
channel  flow  and  transport  phenomena  are  summarized  as  follows: 

1.  Equivalent  ice-covered  and  free-surface  flows  are  considered 
to  be  uniform  two-dimensional  flows  which  have  the  same  unit 
water  discharge,  slope,  and  temperature. 

2.  In  comparison  to  a free-surface  flow,  the  depth  of  an  ice- 
covered  flow  is  increased  due  to  the  drag  on  the  underside 
of  the  ice  cover.  The  amount  of  increase,  which  depends 

on  the  composite  roughness  of  the  ice  cover  and  bed  relative 
to  the  bed  roughness  for  the  equivalent  free-surface  flow,  is 
typically  10  to  40  percent.  There  is  a corresponding  reduc- 
tion in  velocity. 

3.  Flows  with  ice  covers  are  characterized  by  lower  and  upper 
layers  of  respective  thicknesses  Yj  and  Y.,,  with  the  locus  of 
maximum  velocity  and  zero  shear  stress  at  the  interface  be- 
tween the  two  layers.  For  linear  shear  stress  distributions 
in  both  layers,  the  thicknesses  of  the  lower  and  upper  layers 
are  proportional  to  the  relative  magnitudes  of  the  shear 
stresses  at  the  lower  and  upper  boundaries,  and  the  velocity 
distribution  in  both  layers  can  be  represented  by  the  Karman- 
Prandtl  logarithmic  velocity  defect  law.  It  follows  that  the 
velocity  and  shear  stress  distribution  in  the  lower  layer,  and 
possibly  the  sediment -transport  and  bed-form  characteristics  to 
a considerable  extent  also,  should  be  the  same  as  in  a free- 


surface  flow  with  depth  Yj  and  mean  velocity  l^. 


4.  The  two-layer  model  can  also  be  extended  to  estimate  the  ver- 
tical distribution  of  suspended  sediment  In  an  ice-covered 
flow.  It  is  assumed  that  the  vertical  turbulent  mass  transfer 
coefficient  for  suspended  sediment  is  distributed  parabolically 
like  the  eddy  diffusivity  in  the  regions  near  the  lower  and  upper 
boundaries  but  that  it  is  linearly  distributed  between  the  ver- 
tices of  the  two  parabolas  in  the  mid-depth  region. 

5.  From  generalized  sediment  transport  relationships  and  also  from 
elementary  work-energy  considerations,  it  is  seen  that  both 
the  total  bed-material  discharge  and  the  bed-load  discharge 
can  be  greatly  reduced  by  the  presence  of  an  ice  cover.  This 
results  from  the  reduction  in  both  mean  velocity  and  bed  shear 
stress  due  to  the  ice  cover. 

6.  From  a consideration  of  currently  available  relationships  for 
predicting  the  type,  characteristic  dimensions,  and  hydraulic 
roughness  of  alluvial  channel  bed  forms,  it  appears  likely  that 
an  ice  cover  causes  a change  in  bed-form  characteristics  similar 
to  that  which  would  occur  in  a free  surface  flow  if  the  flow 
discharge  were  reduced,  or  if  the  discharge  were  held  constant 
and  the  depth  increased.  For  flows  in  the  ripple  and  dune 
regimes  a reduction  in  bed-form  size,  and  possibly  an  increase 
in  bed  roughness  due  to  decreasing  bed-form  length-to-height 
ratio,  is  expected.  This  is  consistent  with  the  hypothesis 
that  bed  friction  factors,  bed-form  characteristics,  and  sedi- 
ment-transport properties  of  an  ice-covered  alluvial  channel  are 
equivalent  to  those  in  a free-surface  flow  with  depth  Yj  and 


mean  velocity  U^. 


III.  EXPERIMENTAL  INVESTIGATION 


A.  Description  of  Experiments.  The  main  purpose  of  the  ex- 
perimental investigation  was  to  obtain  comparisons  between  the  charac- 
teristics of  alluvial-channel  flows  with  a simulated  ice  cover  and  equiv- 
alent (same  discharge  and  slope)  free-surface  flows.  To  this  end,  three 
sets  of  experiments  were  performed  in  a laboratory  flume  with  a sand 


bed.  Each  set  consisted  of  three  runs,  one  with  a free  surface,  and  one 
each  with  a smooth  and  a rough  simulated  ice  cover.  The  slope  and  dis- 
charge were  supposed  to  be  the  same  for  the  three  runs  in  each  set;  how- 
ever, for  reasons  to  be  explained  later,  equality  of  discharge  was 
achieved  in  only  one  set.  Data  obtained  during  each  run  included  measure- 
ments of  discharge,  water-surface  elevations  and  slope,  water  temperature, 
several  longitudinal  bed  profiles,  and  concentration  of  the  total  sediment 
load.  Velocity  and  suspended  sediment  concentration  profiles  were  ob- 
tained for  all  runs  in  two  of  the  sets.  The  experimental  equipment  and 
procedures  are  described  in  the  next  two  sections. 

1.  Equipment.  The  principal  experimental  facility  was  a re- 
circulating flume  with  a glass-walled  working  section  measuring  27.4  m 
in  length,  0.914  m in  width,  and  0.45  m in  depth.  An  elevation  view  of 
the  flume  is  shown  in  figure  3.  The  power  to  recirculate  the  water-sedi- 
ment mixture  was  provided  by  two  7.5-horsepower  variable  speed  motor-pump 
units  located  under  the  downstream  end  of  the  flume,  each  discharging 
into  separate  10-inch  spiral  weld  return  pipes  terminating  at  a transi- 
tion to  the  flume  inlet.  Discharge  could  be  measured  at  side-contracted 
orifices  in  each  return  pipe.  The  flume  was  supported  by  a central  pivot 
and  four  synchronized  motor-driven  cam  jacks  so  that  the  slope  could  be 
changed  without  interrupting  operation  of  the  flume.  Flume  slope  could 
be  measured  by  means  of  a dial  gage  and  a point  gage  at  the  downstream 
end  of  the  flume  for  measuring  vertical  deflection  from  the  horizontal 
position.  Steel  rails  for  the  instrument  carriage,  mounted  atop  the 
flume  walls,  provideda  reference  frame  for  elevation  relative  to  a plane 
tilted  at  the  flume  slope. 

The  bed  was  composed  of  quartz  foundry  sand  with  a median  diameter, 
D50,  of  0.25  mm  and  a geometric  standard  deviation,  o^,  of  1.41.  A plot 
of  the  sieve-diameter  distribution  for  the  sand  is  shown  in  figure  4. 

The  same  flume  and  sand  were  used  in  experimental  investigations  reported 
by  Onishi  et  al.  (1972),  Jain  (1971),  Annambholta  (1969),  and  Squarer  (1968). 

A smooth  ice  cover  was  simulated  by  1.27-cm  thick,  0.91  m by  1.22  m 
painted  plywood  panels,  connected  by  hinges,  so  that  they  constituted  a 
nonrigid  floating  cover  extending  over  almost  the  entire  length  of  the 
27.4-m-long  working  section.  Each  panel  was  reinforced  to  prevent  warping. 
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To  simulate  a rough  ice  cover,  0.64-cm  thick  by  2.54-cm  wide  continu- 
ous masonite  strips,  extending  across  the  width  of  the  panels,  were 
tacked  to  the  bottom  surface  of  the  panels  at  intervals  of  15.24  cm. 

Water  surface  elevations  and  slopes  were  measured  by  means 
of  piezometers  consisting  of  perforated,  gauze-wrapped  copper  tubes, 
buried  in  the  sand  and  extending  across  the  flume  at  3.048-m  inter- 
vals along  the  working  section  of  the  flume.  The  piezometers  were  con- 
nected by  flexible  plastic  tubing  to  a bank  of  glass  manometer  tubes, 
where  the  water  level  in  each  tube  could  be  read  to  the  nearest  0.3  mm 
by  means  of  a horizontal  wire  attached  to  a point-gage  assembly. 

Longitudinal  bed  profiles  were  measured  with  an  Automation 
Industries  Model  1054  ultrasonic  distance  meter  that  was  mounted  to 
the  motor-driven  instrument  carriage  so  that  the  sensing  face  of  the 
transducer  was  submerged  in  the  water  about  8 cm.  above  the  bed.  The  ana- 
log output  from  the  distance  meter  was  digitized  and  sampled  at  80- 
millisecond  intervals  by  the  Institute's  IBM  1800  Data  Acquisition  and 
Control  System  which  performed  a preliminary  test  for  stationarity  and 
then  punched  the  profile  data  on  IBM  cards  pending  acceptance  of  the 
data  as  satisfactory.  The  carriage  speed  was  set  so  that  the  80-milli- 
second sampling  interval  would  correspond  to  approximately  1-cm  inter- 
vals along  the  working  section  of  the  flume.  The  punched  IBM  cards 
were  subsequently  taken  to  the  IBM  360/65  at  the  University  Computer 
Center  for  further  statistical  analysis  and  plotting  by  an  incremental 
plotter. 

Samples  for  determining  total  sediment  load  concentrations 
were  obtained  by  withdrawing  a representative  water-sediment  mixture 
through  a vertical  slot  sampler,  consisting  of  two  parallel  sheet  metal 
plates  connected  to  a hose  and  pump,  at  the  downstream  end  of  the  flume. 
This  device  sampled  a 1-cm  wide  vertical  slice  extending  over  the  entire 

3 

depth  of  flow  at  the  centerline.  Samples  were  pumped  into  a 0.65-m  col- 
lecting tank  at  such  a rate  that  the  mean  entrance  velocity  into  the 
slot  equalled  the  depth-averaged  flow  velocity  in  the  flume.  Makeup 
water  was  added  during  sampling  to  maintain  a constant  volume  of  water 
in  the  flume  system.  Following  collection  and  sufficient  time  for 
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the  sediment  to  settle,  the  water  was  decanted  and  the  sediment  was 
carefully  flushed  into  a beaker  for  subsequent  drying  and  weighing. 

Velocity  profiles  were  measured  by  means  of  a 
6.3-mm  outside-diameter  Prandtl  Pitot  tube  mounted  to  a point-gage 
assembly  and  connected  to  a dual  pressure-transducer  system.  The  dy- 
namic and  static  pressures  were  transmitted  by  water  columns  to  two  ± 1 kPa 
capacity  Model  PM5-TC  Statham  pressure  transducers.  The  outputs  of  the 
transducers  were  amplified  by  noninverting  high-gain  chopper  amplifiers 
and  fed  sequentially  through  a difference  amplifier  and  voltage-to- 
frequency  converter  to  an  electronic  counter  for  display.  Sensitivity 
balancing  of  the  transducer  outputs  by  the  difference  amplifier  reduced 
errors  due  to  temperature  fluctuations  and  vibrations.  The  system  was 
capable  of  accurately  measuring  velocities  down  to  about  3 cm/sec.  Cali- 
bration was  performed  by  placing  the  Pitot  tube  in  the  center  of  a sub- 
merged jet  just  downstream  from  a flow  nozzle,  the  discharge  of  which 
had  been  determined  by  weighing  the  flow  accumulated  over  a measured 
time  interval. 

Samples  for  determining  suspended  sediment  concentration 
profiles  were  siphoned  through  a sampling  nozzle,  attached  to  a point- 
gage  assembly,  into  400-ml  beakers.  The  intake  end  of  the  sampling 
nozzle  was  shaped  into  a rectangle,  3.2  mm  high  and  11.1  mm  wide.  For 
each  sample,  the  outlet  elevation  of  the  siphon  tube  was  fixed  so  that 
the  entrance  velocity  into  the  sampling  nozzle  would  be  equal  to  the 
local  flow  velocity  which  had  just  been  measured  at  the  sampling 
point.  After  weighing  the  total  samples  and  allowing  the  sediment  to 
settle,  the  excess  water  was  decanted  and  the  sediment  and  remaining 
water  were  transferred  to  evaporation  dishes  for  oven-drying  and  weighing. 
Samples  were  weighed  to  the  nearest  milligram  with  a Sartorius  elec- 
tronic precision  balance. 

2.  Experimental  procedures.  Three  sets  of  three  runs  each 
were  conducted,  constituting  a total  of  nine  runs.  They  were  designated 
by  two  alphabetic  characters:  the  first  (A,B,  or  C)  in  reference  to 
the  slope,  S,  and  unit  discharge,  q,  for  the  set;  and  the  second  (F, 

S,  or  R)  to  denote  flow  with  a free  surface  (F) , a smooth  simulated  ice 


cover  (S),  or  a rough  simulated  ice  cover  (R) . Thus  Run  BR,  for  example, 
means  a run  with  slope  and  discharge  set  B,  and  a rough  simulated  ice 
cover. 

The  original  intent  was  for  the  three  runs  in  each  set  to 
have  equivalent  flow  conditions  in  the  sense  that  they  all  be  character- 
ized by  the  same  slope  and  discharge.  Unfortunately  it  was  discovered 
following  sets  A and  B that  equality  of  discharge  in  these  sets  had 
not  been  achieved.  In  set  C,  the  condition  causing  the  unequal  dis- 
charges was  corrected,  so  that  the  discharges  as  well  as  the  slopes 
were  the  same.  The  difference  in  the  discharges  for  the  three  runs  in 
sets  A and  B occurred  because  only  one  of  the  two  pumps  in  the  parallel 
return-flow  pipe  system  was  operated  during  these  runs  because  it  was 
considered  that  at  low  discharges,  the  discharge  could  be  measured 
more  accurately  and  there  would  be  less  deposition  of  sediment  in  the 
pipe  with  this  mode  of  operation.  However,  since  there  are  no  shutoff 
valves  in  the  system  - discharge  control  is  achieved  entirely  by  regu- 
lating pump  speed  - a portion  of  the  flow  could  be  recirculated  from 
the  headbox  of  the  flume  back  through  the  other  pipe  instead  of  through 
the  flume.  To  the  extent  that  this  occurred,  the  discharge  in  the  flume 
was  less  than  that  indicated  by  the  orifice  and  manometer  unit  for 
the  pipe  with  the  operating  pump.  Comparison  of  discharges  computed 
by  integrating  the  measured  velocity  profiles  over  the  flow  cross 
section  in  the  flume  and  those  indicated  by  the  orifice  and  manometer 
unit  showed  that  the  backflow  through  the  second  pipe  was  negligible 
in  the  runs  with  a free  surface,  but  significant  in  the  runs  with  a 
simulated  ice  cover  due  to  the  added  resistance  to  flow  through  the 
flume.  Consequently  the  discharges  reported  for  sets  A and  B are  the 
ones  which  were  computed  by  integrating  the  measured  velocities  over  the 
flow  cross  section. 

The  experimental  procedures  were  somewhat  different  for  the 
free-surface  runs  and  the  runs  with  the  simulated  ice  covers.  The  free- 
surface  run  was  always  the  first  of  a set.  Before  starting  the  run, 
the  sand  bed  was  leveled  over  the  length  and  width  of  the  flume  with 
the  aid  of  a screed  that  was  mounted  to  the  instrument  carriage.  Next, 
with  the  flume  slope  horizontal,  enough  water  was  added  to  the  flume  to 


bring  the  depth  to  an  estimated  0.100  m.  Then  the  flume  slope  was  set 
at  the  estimated  value  for  uniform  flow,  and  the  flow  was  started  at 
the  desired  discharge.  From  time  to  time  the  water  surface  slope  was 
checked  by  reading  the  water-surface  elevations  on  the  manometer  bank, 
and  the  flume  slope  was  adjusted  as  necessary  to  make  it  equal  to  the 
water-surface  slope.  Discharge  would  also  be  checked  and  any  necessary 
adjustments  made  at  these  times.  This  procedure  of  successive  adjust- 
ments was  continued  as  the  bed-form  configuration  came  into  equilibrium 
with  the  flow  conditions.  After  the  water-surface  slope  had  remained 
equal  to  the  flume  slope  for  a period  of  several  hours,  usually  over- 
night, and  the  flow  depths  at  several  points  along  the  glass  walls  had 
been  checked  with  a tape  measure  to  see  if  the  flow  depth  was  approxi- 
mately 0.10  m along  the  entire  length  of  the  flume,  the  flow  was  stopped 
and  the  bed  profile  along  the  centerline  was  measured  with  the  ultra- 
sonic distance  meter.  The  output  from  the  IBM  1800  Data  Acquisition 
and  Control  System  was  reviewed  for  nonstationarity , i.e.  any  consistently 
rising  or  falling  trend  along  the  flume  in  the  bed  elevation  with  respect 
to  the  mean  bed  elevation.  If  a visual  scan  indicated  stationarlty  of 
the  data,  the  slope,  discharge  and  average  depth  were  recorded  and  the 
bed  profile  data,  stored  on  disk  by  the  computer,  were  punched  on  IBM 
cards;  otherwise  the  flow  was  started  again,  and  the  process  was  repeated 
after  a few  more  hours.  If  the  centerline  data  passed  the  stationarity 
test,  bed  profiles  were  obtained  along  two  additional  sailing  lines, 
one  quarter  of  the  flume  width  to  the  left  of  and  to  the  right  of  the 
centerline.  The  flow  was  then  started  again  and  five  or  six  total  sedi- 
ment load  samples  were  obtained,  taking  care  to  add  sufficient  makeup 
water  to  maintain  uniform  flow  at  the  proper  depth  and  slope.  Then 
three  sets  of  velocity  and  suspended  sediment  concentration  profiles 
were  obtained  at  about  mid-length  along  the  working  section,  one  each 
at  the  centerline  and  one-quarter  flume  width  to  the  left  and  right 
of  the  centerline.  Velocity  measurements  were  obtained  at  twelve  to 
fifteen  points  in  the  vertical,  and  suspended  sediment  samples  at  five 
points.  The  flow  was  then  stopped  again  and  a second  set  of  three  longi- 
tudinal bed  profiles  was  obtained,  along  the  centerline  and  along  the 
left  and  right  sailing  lines.  Water  temperature  was  measured  with  a 
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mercury  thermometer  at  intervals  throughout  the  run.  This  completed 
the  procedure  for  a free-surface  run. 

Except  for  the  differences  noted  below,  the  procedures  for 
the  runs  with  the  simulated  ice  covers  were  the  same.  In  these  runs, 
the  flume  slope  was  maintained  constant  at  the  value  determined  for  the 
corresponding  free  surface  run,  i.e.  the  first  run  of  the  set.  The  bed 
was  not  leveled  and  screeded  before  installing  the  simulated  ice  cover 
and  starting  the  flow  at  the  beginning  of  the  run.  During  the  first 
part  of  the  run,  the  amount  of  water  in  the  flume  was  increased  and 
adjusted  by  successive  trials  until  the  piezometric  slope  became  and 
remained  equal  to  the  flume  slope.  Before  obtaining  the  longitudinal 
bed  profiles  and  after  stopping  the  flow,  the  panels  simulating  the  ice 
cover  had  to  be  removed.  This  was  done  with  care  to  minimize  disturbance 
to  the  bed  forms.  The  panels  were  replaced  before  obtaining  the  total 
sediment  load  samples  and  velocity  and  the  suspended  sediment  concentra- 
tion profiles.  Then  the  panels  had  to  be  removed  once  more  before  obtaining 
the  final  set  of  bed  profiles. 

In  set  C the  procedures  were  varied  somewhat  from  those  outlined 
above.  Both  pumps  in  the  flume  recirculating  system  were  operated 
which  eliminated  the  problem  of  backflow  through  the  second  pipe.  Also, 
no  velocity  and  suspended-sediment  concentration  profiles  were  obtained, 
and  bed  profiles  were  taken  only  along  the  centerline. 

B.  Summary  of  Experimental  Results.  The  bulk  flow  and 
sediment-transport  characteristics  for  each  run  are  listed  in  table  1. 

All  quantities  listed  are  averages  obtained  in  one  or  more  of  the  following 
ways:  averages  of  measurements  taken  at  various  times  throughout  the 
runs,  averages  of  measurements  obtained  at  the  centerline  and  at  the 
left  and  right  quarter  points  across  the  flume,  averages  along  the  length 
of  the  working  section,  and  averages  over  the  depth  of  flow.  The  unit 
discharge,  q,  was  determined  from  the  average  of  the  centerline  and  the 
left  and  right  quarter-point  velocity  profiles  for  the  series  A and  B 
runs,  and  from  the  orifice-meter  discharge  measurements  for  the  series 
C runs.  The  mean  depth,  d,  was  determined  by  subtracting  the  mean  bed 
elevation  from  the  mean  water-surface  elevation,  where  the  former  was 
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determined  from  the  bed  profiles  along  the  centerline  and  the  left 
and  right  sailing  lines,  and  the  latter  from  the  piezometer  readings 
along  the  working  section  of  the  flume.  The  mean  velocity,  0,  is  equal 
to  q/d.  The  slope,  S,  is  the  flume  slope  which  in  turn  was  set  to 
equal  the  water-surface  slopes  indicated  by  the  piezometers.  The  tempera- 
ture, T,  and  mean  total  sediment  concentration,  C^,  are  averages  of 
measurements  and  samples  taken  several  times  during  the  course  of  the 

runs.  The  unit  total  sediment  discharge,  q is  equal  to  qC  multiplied 

kg  l St  -3  C 

by  the  appropriate  conversion  factor  (— “ — ~ =10  ) . 

m3 

The  bed  configuration  consisted  of  ripples  for  all  runs.  A 
statistical  analysis  of  the  bed  forms  is  presented  in  section  111-h. 

Data  showing  the  vertical  distributions  of  velocity,  u,  sus- 
pended sediment  concentration,  C,  and  the  suspended  sediment  flux  uC, 
for  the  A and  B series  of  runs,  are  presented  in  figure  5,  (a)  through 
(f).  The  curves,  which  are  based  on  the  average  of  the  three  measure- 
ments taken  at  the  centerline  and  the  right  and  left  quarter  points,  were 
drawn  with  some  smoothing. 

Various  quantities  determined  from  the  velocity  profiles  are 
listed  in  table  2.  The  reference  depths  and  velocities,  d and  U,  for 
the  profiles  in  figure  5,  are  averages  for  the  three  verticals  at  which 
the  profiles  were  measured.  The  values  of  d and  U,  which  are  listed 
in  table  2,  are  in  general  not  equal  to  the  mean  values,  d and  0,  for  the 
entire  working  section  listed  in  tabic  1. 

For  both  the  velocity  and  concentration  data  in  figure  5, 
departures  from  the  mean  profile  appear  to  be  more  systematic  than  random. 

For  the  velocity  data,  the  tendency  for  departures  to  be  greater  near 
the  bed  is  likely  due  to  variations  in  the  positions  of  the  Pitot  tube 
relative  to  the  bed  forms.  For  the  suspended  sediment  concentration 
data,  there  is  a fairly  consistent  tendency  for  the  concentrations 
measured  at  the  left  and  right  one-quarter  points  to  be  fairly  close  to 
each  other,  but  for  the  concentrations  measured  at  the  centerline  to  be 
either  significantly  higher  or  lower  - more  often  higher  in  the  free- 
surface  runs,  and  lower  in  the  runs  with  the  simulated  ice  cover.  This 
could  be  due  to  the  existence  of  secondary  circulation  cells,  wherein 
an  upward  velocity  component  would  tend  to  increase  the  local  concentration 


uC/U  in  mg/L 


(c)  Run  AR 

Figure  5.  (continued). 


ative  Depth, 


Table  2.  FLOW  CHARACTERISTICS  DETERMINED  FROM  VELOCITY  MEAS 
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and  a downward  component  to  decrease  it.  This  tendency  is  not  evident 
in  the  velocity  data.  In  keeping  with  the  above  observations,  it  is 
to  be  expected  that  velocities  would  be  relatively  more  responsive  to 
local  stimuli  such  as  the  location  of  a bed  form  with  respect  to  the 
Pitot  tube,  whereas  the  local  concentrations  would  depend  to  a greater 
extent  on  the  integrated  effects  of  mass  transfer  phenomena  which  have 
occurred  for  some  distance  upstream. 

The  mean  suspended  sediment  concentrations,  Cs>  listed  in 
column  10  of  table  1 were  computed  by  integrating  the  average  sediment- 
flux  distribution  (uC/U  vs  q)  curves  in  figure  5,  over  the  depth  of 
flow.  The  unit  suspended  load  discharge,  qgs»  is  equal  to  qC^  multiplied 
by  the  appropriate  conversion  factor.  No  satisfactory  explanation  is  at 

hand  as  to  why  the  C values  exceed  the  C values  in  Runs  AS,  AR  and  BF. 

J s t 

Overall,  the  C values  are  considered  to  be  more  reliable  than  the  C 
t s 

values.  However,  the  variability  of  individual  values  from  the  average 
values  listed  in  table  1 for  the  A and  B series  of  runs  (ranging  up  to 
± 50%  and  more)  was  much  greater  than  the  variability  for  the  C series. 

This  suggests  that  an  unstable  condition  due  to  the  recirculation  back 
through  the  second  return  flow  pipe  may  have  existed  at  the  downstream 
end  of  the  flume  near  the  slot  sampler  and  affected  the  total-sediment 
concentration  samples. 

The  quantities  listed  in  columns  5 through  13  of  table  2 
were  determined  in  the  process  of  fitting  (6)  to  the  velocity-distribution 
data  in  figure  5.  They  will  be  discussed  in  the  next  section. 

Average  velocities,  shear  velocities,  and  Manning  and  Darcy- 
Weisbach  friction  factors  for  the  composite  channel,  and  for  the  lower 
and  upper  layers  are  listed  in  table  3.  They  were  computed  from  the  bulk 
flow  characteristics  in  table  1,  using  the  Y^/d,  Y^/d,  Uj/U,  and  U2/U  values 
in  table  2 that  were  determined  from  the  velocity  measurements.  It 
is  assumed  that  these  ratios  remained  constant  along  the  length  of  the 
working  section.  Values  of  the  composite  shear  velocity  were  computed 
from  the  equation 


Table  3.  AVERAGE  VELOCITIES,  SHEAR  VELOCITIES,  AND  FRICTION  FACTORS 
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and  values  of  n, , f,t  u.  , and  u.  , n,  and  n„,  and  f.  and  f„,  according 
1 I *2  l ^ i z 

to  (17),  (.18),  (8),  (15),  and  (16)  respectively.  These  results  will  be 
discussed  in  the  sections  on  bed  friction  factors  and  depth  of  flow. 


C.  Velocity  Distribution.  Using  the  relative  depth,  n = y/d, 
where  y is  measured  from  the  bed,  instead  of  y^  and  y^  as  defined  in 
figure  1,  as  independent  variables,  the  velocity  distribution  according 
to  (6)  becomes 


0 < r)  < Y1/d: 


Y /d  £ n £ 1: 


u.  - u 
1 max 


ln[n  d/Y  ] 

*1  1 


ln[(l-n)d/Y2] 


Figure  6 shows  (46)  fitted  to  the  average  velocity  profile  data  for 

each  run.  The  fitting  was  accomplished  in  the  following  manner. 

After  first  averaging  the  velocities  measured  at  corresponding  distances 

below  the  upper  surface  in  the  verticals  at  the  centerline  and  the  left 

and  right  quarter  points,  the  average  velocities  were  plotted  as  a 

function  of  n = y/d,  and  smooth  curves  were  drawn  by  eye  through  the 

plotted  points.  In  the  averaging,  velocities  for  points  below  the 

local  bed  level  were  assigned  a value  of  zero.  Values  of  Y,/d  and  u 

were  then  read  from  the  curves,  and  values  of  U,  U^,  and  were  computed 

by  numerically  integrating  the  curves.  Values  of  u.  were  then  computed 

1 

by  (8),  with  Y^  = d (Y^/d)  and  d as  the  average  depth  for  the  three 
verticals.  Finally,  values  of  and  were  computed  by  (14),  completing 
the  information  needed  for  computing  and  plotting  the  curves  in  figure  6. 

The  values  of  the  various  parameters  determined  during  the 
fitting  process  are  listed  in  columns  5 through  13  of  table  2.  These 
results  show:  (1)  that  Y^/d  appears  to  depend  primarily  on  the  roughness 
of  the  ice  cover,  and  to  be  relatively  insensitive  to  q,  implying  that 
the  bed  roughness,  given  a particular  ice  cover,  may  also  be  relatively 
insensitive  to  q;  (2)  that  the  ratio  u /U  remains  relatively  constant 
for  different  combinations  of  upper  and  lower  boundary  roughnesses  and 
flow  conditions;  and  (3)  that  the  mean  velocities  in  the  lower  and  upper 
layers,  U^  and  U2>  are  quite  close  to  the  overall  mean  velocity  0.  In 


RUN  AF  # RUN  AS  f RUN  AR 


showing  (46)  fitted  to  velocity  distribution  data 


42 


addition,  there  is  considerable  variation  in  the  k values;  in  particular 

those  for  the  upper  layer  with  the  rough  cover  seem  to  be  very  large. 

However,  the  < values  computed  by  (13)  were  found  to  be  in  reasonably 

good  agreement  with  ones  determined  more  conventionally  from  semilog 

plots  of  u/u^  vs  y /Y . , although  the  latter  tended  to  be  slightly 
i 11 

lower. 


D.  Depth  of  Flow.  The  results  in  column  3 of  table  1 
clearly  show  that  the  mean  depth  of  flow  is  increased  by  an  ice  cover. 
However,  because  of  the  reduction  in  discharge,  due  to  recirculation 
back  through  the  second  return-flow  pipe,  which  occurred  in  the  series 
A and  B runs  with  simulated  ice  covers , the  depth  increases  observed  in 
these  runs  were  not  the  same  as  they  would  have  been  had  the  discharges 
been  the  same  as  in  the  corresponding  free-surface  runs.  The  following 
analysis  of  the  variation  of  depth  with  discharge  in  an  ice-covered 
alluvial  channel  provides  a method  for  estimating  the  depths,  which 
would  have  occurred  for  equivalent  ice-covered  flows  in  Runs  AS,  AR,  BS, 


and  BR. 


The  relationship  between  the  depth  and  the  other  bulk  flow 


characteristics  for  an  ice-covered  channel  is 


_.-l/3  2/3  . 1/3 

(4gS)  q f 


according  to  the  Darcy-Weisbach  formula,  (18),  and 

. 90.4  -0.3  0.6  0.6  ,/Q. 

dj.  = 2 S q n-j.  (48) 

according  to  the  Manning  formula,  (17).  In  response  to  a change  in  q 
in  a sand-bed  channel,  not  only  d^,  but  also  the  bed  friction  factors, 
f^  and  n^,  and  consequently  the  composite  friction  factors,  f^  and  n^ , 
are  subject  to  change.  Approximate  expressions  for  the  composite  friction 
factors  in  terms  of  the  bed-friction  factors  and  Y^/d  result  from  com- 
bining (20)  and  (26)  to  obtain 


2Y1/d 


and  (19)  and  (25)  Co  obtain 


Hj.  = ^/d)”273  nx 


(50) 


Finally,  in  keeping  with  the  hypothesis  that  flow  in  the  lower  layer 
resembles  a free-surface  flow  with  depth  and  mean  velocity  U^,  the 
experimentally-determined  values  of  f^  and  n^  are  plotted  in  figure  7 
as  functions  of  unit  discharge  in  the  lower  layer 

'll  ‘ Vi  * T ’ <51) 

It  should  be  noted  that  the  experimental  conditions  for  the  data  in 
figure  7 are  restricted  to  a narrow  range  of  slopes  (0.0018  to  0.0021) 
and  a single  set  of  bed  material  characteristics . In  using  (49),  (50), 
and  (51)  it  is  also  implied  that  = U,  and  that  Y^/d  is  inde- 

pendent of  q.  The  data  in  table  2 appear  to  indicate  that  both  of  these 
conditions  are  satisfied,  at  least  approximately. 

Estimated  values  of  the  depth,  d /d  , for  ice-covered  flows 

X F 

that  are  equivalent  to  the  free-surface  flows  in  Runs  AF,  BF,  and  CF 
are  listed  in  table  4.  The  values  of  3^.  in  columns  7 and  9 were 
computed  by  (47)  and  (48) , using  the  q and  S values  in  columns  2 
and  3 , and  the  f^  and  n^  values  from  figure  7 in  (49)  and  (50)  to 
obtain  the  f^.  and  n^.  values.  The  small  differences  between  the  esti- 
mates based  on  the  Darcy  Weisbach  and  Manning  formulas  are  not  consid- 
ered to  be  significant.  The  d^  values  in  column  (5)  represent  the 
experimentally  observed  depths.  They  should  be  compared  directly 
with  the  dj.  estimates  in  columns  7 and  9 only  for  the  C series  of 
runs  wherein  the  discharge  remained  constant.  Such  comparison  indi- 
cates quite  good  agreement  between  the  observed  d values  and  the  d 

A 

estimates  for  Runs  CS  and  CR.  Differences  between  d^  and  d^  for  the  A 
and  B series  of  runs  are  of  course  due  mainly  to  the  recirculation- 
related  discharge  reductions  in  the  experiments. 


E.  Bed  Friction  Factors.  The  hypothesis  that  the  flow  in  the 


lower  layer  of  an  ice-covered  flow  is  like  that  in  a free-surface  flow  of 
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depth  Y^  and  mean  velocity  is  tested  further,  using  modified  forms 
of  the  Alam  and  Kennedy  (1969)  bed-form  friction-factor  chart. 

Following  Taylor  and  Brooks  (1961) , the  Darcy-Weisbach  fric- 
tion factor  for  the  bed  is  separated  into  two  parts, 

f1  = + f'l  (52) 

where  fj  is  due  to  the  roughness  of  the  sand  grains,  and  f^  is  due 

to  the  form  drag  associated  with  the  large-scale  features  of  the  bed 

forms.  Values  of  the  bed-form  friction  factor,  f!^,  are  determined  by 

subtracting  f|  from  f ^ , where  fj  is  evaluated  by  means  of  the  Lovera  and 

Kennedy  (1969)  flat  bed  friction  factor  diagram,  where  f j is  represented 

as  a function  of  the  Reynolds'  number,  ]R  = Ur^/v,  and  the  relative 

roughness,  r^/D  q.  The  Alam  and  Kennedy  bed-form  friction-factor  chart, 

which  is  based  on  river  and  flume  data  from  many  sources,  consists  of 

a set  of  empirical  curves  representing  constant  values  of  U//gD^g  on  a 

graph  of  f^  vs  r^/D^g.  For  the  purposes  of  the  present  investigation 

when  using  these  charts,  r.  , the  hydraulic  radius  of  the  bed  is  re- 

b 

placed  by  for  flows  with  simulated  ice  covers,  and  Y^  = d for  free- 

surface  flows.  It  is  furthermore  assumed  that  = U.  Also,  no  sidewall 

corrections  are  made  insofar  as  r,  /d  > 0.90  for  all  free-surface  data 

b 

included  in  the  present  analysis,  and  the  need  for  corrections  would 
be  even  less  for  data  from  experiments  with  simulated  ice  covers.  Values 
of  the  flat-bed  and  bed-form  friction  factors  and  related  quantities  for 
the  present  experiments  are  listed  in  table  5.  The  rows  with  a run 
designation  followed  by  a prime,  e.g.  AS',  contain  values  of  the  various 
quantities  computed  for  equivalent  ice-covered  flows  with  the  d^  values 
in  table  4. 

Figures  8 and  9 show  data  from  the  present  investigation,  to- 
gether with  data  from  Annambhotla  (1969)  and  Onishi  et  al.  (1972)  for 
free-surface  flows  in  the  same  flume  with  the  same  bed  material,  super- 
imposed on  modified  forms  of  the  Alam  and  Kennedy  bed-form  friction-factor 
chart.  The  numbers  next  to  the  plotted  points  denote  values  of  the 
third  variable:  0/ /gD,_g  in  figure  8,  and  in  figure  9.  The  only 


Table  5.  FLAT-BED,  AND  BED-FORM  FRICTION  FACTORS 
AND  RELATED  QUANTITIES 


Run 

(1) 

fl 

(2) 

* - 5Yl 

VD50 

(4) 

fi 

(5) 

U 

fI 

(7) 

V 

(3) 

/gD5o 

(6) 

AF 

0.106 

36,300 

424 

0.018 

7.58 

0.088 

AS 

0.149 

26,100 

378 

0.020 

6.09 

0.129 

AS’* 

0.149 

29,300 

408 

0.018 

6.32 

0.131 

AR 

0.108 

15,000 

218 

0.025 

5.54 

0.083 

AR' 

0.119 

18,000 

454 

0.024 

5.70 

0.095 

BF 

0.057 

54,000 

408 

0.020 

10.93 

0.037 

BS 

0.159 

30  ,900 

390 

0.019 

6.48 

0.140 

BS’ 

0.099 

43,900 

424 

0.019 

8.44 

0.080 

BR 

0.112 

16,300 

237 

0.024 

5.83 

0.088 

BR’ 

0.124 

21,000 

296 

0.023 

6.02 

0.101 

CF 

0.077 

43,600 

428 

0.018 

9.37 

0.059 

CS 

0.113 

35,500 

424 

0.018 

7.52 

0.095 

CR 

0.105 

17,700 

259 

0.024 

6.20 

0.081 

Primes  following  run  designation  identify  rows  corresponding  to 
computed  equivalent  ice-covered  flows  based  on  d values  in 
table  4. 
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Figure  8.-  Modified  Alam  and  Kennedy  bed-form  friction  factor  chart;  f"  vs, 
Y /D  with  u//gD  as  third  variable. 


Figure  9.-  Modified  Alam  and  Kennedy  bed-form  friction  factor  chart 
U//rDTT  with  Y,/D^_  as  third  variable. 
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substantive  modifications  in  figures  8 and  9,  in  comparison  to  the  ori- 
ginal  Alam  and  Kennedy  chart,  are:  (1)  the  substitution  of  Y^  or  d 
for  r,  ; (2)  the  interchanging  of  Y. /D  with  0/ /gD  as  the  abscissa 

D 1 jU  3U 

and  third  variables  in  figure  9;  (3)  restriction  to  a range  of  Y^/D^^ 
and  U//gDj-Q  values  that  is  sufficient  to  include  the  data  to  be  con- 
sidered; and  (4)  representation  in  linear  instead  of  logarithmic  coordi- 
nates. No  numerical  adjustments  were  made  in  transcribing  the  curves 
from  the  original  chart  to  the  formats  shown  in  figures  8 and  9. 

The  data  points  from  the  present  investigation  agree  with 
the  curves  about  as  well  as  the  data  of  Onishi  et  al. , and  the  data 
included  by  Alam  and  Kennedy  on  the  original  chart,  and  better  than  the 
data  of  Annambhotia.  The  points  representing  the  computed  equivalent 
ice-covered  flows,  indicated  by  a prime  following  the  run  designation, 
agree  as  well  with  the  curves  as  do  the  points  representing  the  measured 
data  for  the  same  runs.  Comparison  between  figures  8 and  9 and  figure 

7 shows  that  the  data  from  the  present  investigation  are  arranged  in  a 
similar  pattern  in  all  three  figures.  The  reason  becomes  obvious 
when  it  is  noted  that  the  product  of  the  abscissa  variables  in  figures 

8 and  9 is  equal  to  q^//g  , which  if  plotted  against  f^',  would 

constitute  a dimensionless  version  of  figure  7.  All  in  all,  the  results 
shown  in  figures  8 and  9 are  consistent  with  the  hypothesis  that  Y^ 
bears  the  same  relationship  to  the  bed-form  friction  factor  in  an 
ice-covered  alluvial  channel  flow  as  does  the  full  depth,  d,  in  a free- 
surface  flow. 

F.  Suspended  Sediment.  In  order  to  compute  the  suspended 
sediment  concentration  profiles  given  by  (31),  (33),  and  (35),  so  that 
they  can  be  compared  with  the  measured  profiles  in  figure  5,  it  is  first 
necessary  to  evaluate  the  coefficients  and  8^  from  the  data.  Values 
of  z^  were  obtained  by  determining  the  average  slope  of  log-log  plots 
of  C vs  (Y^  - Yi)/yi»  where  C was  taken  from  the  curves  in  figure  5, 
and  y^  is  as  defined  in  figure  1.  Such  a plot,  for  Run  BS,  is  shown  in 
figure  10(a).  Insofar  as  z^  and  z^  represent  the  theoretical  slopes  of 
the  log-log  plots  only  in  the  bottommost  and  topmost  regions  within  a 
distance  of  Y^/2d  from  the  boundaries  in  figure  2 where  (31)  and  (35) 
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apply,  the  upper  parts  of  the  curves  in  figure  10(a),  where  (Y^  - y^)/y^  > 1» 
should  be  weighted  more  heavily  than  the  lower  parts  in  determining  z^ 
values  for  flows  with  ice  covers.  For  free-surface  flows,  wherein  (31) 
represents  the  concentration  profile  over  the  entire  depth  of  flow, 
the  lower  and  upper  parts  of  the  curve  should  be  weighted  equally. 

Having  determined  z^  and  z^.  8^  and  8^  “ere  then  determined  from  the 
definition  equation  for  z^. 


z . 
i 


eiKi  U*i 


(53) 


given  immediately  following  (31)  and  (35),  using  k and  u.  values 
from  table  2,  and  w^  = 0.033  m/sec  (Interagency  Committee  1957)  for 
naturally  worn  quartz  particles  with  a shape  factor  of  0.7,  a sieve 
diameter  of  0.25  mm,  and  a temperature  of  20°C.  Values  of  z^  and  8^ 
determined  from  the  experimental  data  are  listed  in  table  6. 

Most  of  the  8^  values  listed  in  table  6 appear  to  be  unrea- 
sonably large.  For  the  most  part  they  tend  to  be  larger  for  the  A 

than  for  the  B series,  and  for  the  upper  than  for  the  lower  layer, 

situations  wherein  the  quantity  of  suspended  sediment  in  the  samples 
tended  to  be  smaller  and  probably  less  accurately  determined.  Use  of 
fall  velocities  corresponding  to  the  median  diameter  of  the  bed  material, 
which  may  be  significantly  larger  than  the  fall  velocities  of  the  par- 
ticles that  were  actually  in  suspension,  could  also  cause  $ to  be 
overestimated  by  (53).  Unfortunately,  no  particle-size  distribution 
analyses  were  performed  on  the  suspended  sediment  samples. 

The  proposed  model  for  the  vertical  distribution  of  the  tur- 
bulent mass  transfer  coefficient  for  sediment  illustrated  in  figure 
2 and  formulated  in  (30),  (32),  and  (34),  was  tested  by  comparing 

normalized  distributions  of  e , determined  from  the  data,  with  theoreti- 

s 

cal  distributions  of  the  eddy  diffusivity  e . Normalized  distributions 
of  eg  were  estimated  from  data-based  plots  of  log  C versus  q,  as 
shown  in  figure  10(b)  for  Run  BS,  together  with  (27)  in  the  form 


8,e 

i mo 


w d 
-P. 

8 . e 
i mo 


C_ 

dC 

dn 


w d 

e Id 

8.e  dlnC 
i mo 


(54) 


Table  6.  QUANTITIES  ASSOCIATED  WITH  SUSPENDED 
SEDIMENT  CONCENTRATION  PROFILES 


Run 

(1) 

i 

(2) 

z 

i 

(3) 

6. 

l 

(4) 

G 

mo 

2, 

m /sec 
(5) 

e 1 

es  2 

6.  1 
l mo 

(6) 

^mo 

(7) 

AF 

i 

0.47 

3.6 

0.000322 

1.54 

3.60 

AS 

i 

0.23 

5.6 

2 

0.31 

14 

0.000306 

1.09 

5.92 

AR 

1 

'0.28 

6.6 

2 

0.14 

8.3 

0.000247 

1.23 

9.51 

BF 

1 

0.90 

1.7 

0.000412 

1.17 

1.70 

BS 

1 

0.46 

3.4 

2 

0.13 

25 

0.000262 

1.29 

4.96 

BR 

1 

0.86 

1.9 

2 

0.60 

2.2 

0.000299 

1.32 

2.54 

from  figure  11 
2 from  (36)  and  (56) 


approximated  by 


w d 


8 . e 

i mo 


- 0.4343 


An, 


8.c  1 log  cycle 

i mo  ° 

of  C 


(53) 


In  (55)  e is  the  depth-averaged  eddy  diffusivity,  obtained  by 
mo 

averaging  (13)  over  d = to  get 


1,2 


2,2, 


f I<!u*  + K2U*2(6^)  1 


(56) 


and  An^  ^ ^ cycle  t*le  c^an8e  in  n corresponding  to  the  change  in 

of  C 

one  log  cycle  of  C of  the  tangent  to  the  log  C versus  n curve  at  the 

value  of  n for  which  the  local  value  of  e /fl  e is  being  determined. 

s i mo  ° 

Vertical  distributions  of  e /B.e  estimated  in  the  above 

s i mo 

manner  for  each  run  are  shown  as  the  solid  curves  in  figure  11,  where 

they  are  compared  with  distributions  of  e^/emo»  computed  according  to 

(13),  represented  by  the  dashed  curves.  The  results  in  figure  11  are 

somewhat  inconclusive,  mainly  because  the  local  values  of  e^  obtained 

by  graphical  differentiation  are  so  sensitive  to  small  differences  in 

the  way  the  curves  like  those  in  figure  10(b)  are  drawn.  Considering 

this,  the  c /B.e  distribution  curves  for  Runs  AS,  AR,  BF  and  BS , are 
s i mo 

reasonably  consistent  with  the  proposed  model.  Values  of  e /6,e  , which 

s i mo 

is  equal  to  the  ratio  of  the  area  under  the  solid  curves  to  the  area 
under  the  dashed  curves,  are  listed  in  column  6 of  table  6.  Values 
of  es/£mo  computed  from  (36)  and  (56)  are  listed  in  column  7 . 

Figure  12  shows  suspended  sediment  concentration  profiles 
computed  by  (31),  (33),  and  (35),  compared  with  the  average  of  the  con- 
centrations measured  at  the  centerline  and  the  right  and  left  quarter 
points.  The  agreement  in  all  cases  is  quite  good.  It  follows  from  the 
extreme  sensivity  of  the  distribution  of  cg  to  small  changes  in  concen- 
tration profiles,  referred  to  in  connection  with  figure  11,  that  computed 
concentration  profiles  are  relatively  insensitive  to  quite  large 
variations  in  the  distribution,  if  not  the  magnitude,  of  e^. 
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For  all  of  the  computed  profiles  in  figure  12,  the  concentra- 
tion at  n = Y /2d  indicated  by  the  corresponding  measured  profile  was 

cl  i- 

used  as  the  reference  concentration,  C . Efforts  to  predict  reference 

SL 

concentrations  for  conditions  characterizing  the  experiments  were  un- 
successful. Predictions  based  on  the  level-of-maximum-f lux  method 
suggested  by  Willis  and  Kennedy  (1977)  led  to  reference  concentrations 
that  are  too  high  by  about  two  to  three  orders  of  magnitude.  Estimates 
based  on  predictions  of  the  bed-load  discharge  according  to  the  modified 
Einstein  procedure  (Colby  and  Hembree,  1955)  resulted  in  reference 
concentrations  that  are  from  about  1.5  to  30  times  too  high. 


G.  Total  Bed-Material  Discharge.  From  the  unit  total 
sediment  discharge,  q , values  listed  in  table  1,  it  is  clear  that  the 
sediment  discharge  in  flows  with  simulated  ice  covers  was  always  con- 
siderably less  than  in  the  corresponding  free  surface  flows.  However, 
table  1 does  not  indicate  how  much  of  the  reduction  is  due  to  the  cover 
and  how  much  is  related  to  the  decrease  in  flow  discharge  due  to  back- 
flow  through  the  second  return  pipe  in  the  A and  B series  experiments. 
The  following  analysis  was  undertaken  with  a view  toward  isolating  ice- 
cover  from  the  other  causative  factors. 

Figure  13  is  based  on  an  empirical  relationship  obtained  by 
Onishi  et  al.  (1972)  from  dimensional  analysis  and  data  for  free-surface 
flows  in  the  same  flume  with  the  same  sand  that  was  used  in  the  present 
investigation.  The  format  of  figure  13  differs  from  the  original  in  the 
substitution  of  for  r^  and  the  adoption  of  0.18  m/s ec  instead  of  0.24 
m/sec  for  U^,  the  critical  velocity  for  incipient  motion.  The  lower 
value  for  is  more  consistent  with  the  results  compiled  by  Hjulstrom 
(1935)  and  the  observation  of  several  investigators  that  the  critical 
velocity  for  a rippled  bed  is  substantially  less  than  for  a flat  bed. 

The  equation  for  the  curve  is 


st 


»u\w 


U-U  . 
0.03  [ ~]A 


(57) 


If  d were  used  instead  of  Y^,  the  points  representing  runs  with  simu- 
lated ice  covers  would  be  displaced  downward  and  to  the  left  almost 


— 


j 


i 


parallel  to  the  curve,  so  it  makes  little  difference  whether  Y^  or  d is 

used.  As  in  the  original  plot  presented  by  Onishi  et  al.,  inclusion  of 

Yl^°50  °r  ^^D50  3S  3 variable  does  not  improve  the  relationship 

perceptibly.  Evidently  the  total  sediment  discharge  depends  mainly  on 

velocity  and  is  relatively  insensitive  to  depth,  at  least  over  the  limited 

range  of  depths  included  in  the  two  sets  of  data  in  figure  13. 

The  primary  dependence  of  total  sediment  discharge  on  velocity 

is  confirmed  by  figure  14,  wherein  q is  plotted  against  U-U  alone, 

St  c 

with  hardly  any  increase  in  scatter  in  comparison  to  figure  13.  The 
equation  of  the  curve  in  figure  14  is 

qst  = 2<u-u/  (58) 

In  figure  15  q is  represented  as  a function  of  Ut^,  the 
unit  stream  power  exerted  on  the  bed,  in  accordance  with  the  simple 
work-energy  analysis  presented  in  (41)  through  (43) . The  equation  of  the 
curve  in  figure  15  is 

qst  = 0.02[Ut0  ]4.  (59) 

The  amount  of  scatter  in  figure  15  is  not  significantly  more  than  in 

figures  13  or  14.  However,  a review  of  the  data  does  not  indicate  that 

there  is  any  consistent  dependency  on  slope  independently  of  U and  Y^. 

The  relationships  in  figures  13,  14,  and  15  were  used  to  obtain 

the  estimates  of  reduction  in  total  sediment  discharge  due  to  ice  cover 

that  are  shown  in  columns  5,  6,  and  7 of  table  7.  These  estimates  were 

made  by  computing  the  q ^/q  rat*os  according  to  (57)  (58),  and  (59), 

using  the  d and  d values  in  table  4 to  compute  the  U and  Y.  values 

needed  in  the  equations.  The  reductions  computed  according  to  (57)  and 

(58)  agree  fairly  closely  with  one  another,  but  not  very  well  with  those 

computed  according  to  (59).  The  values  of  q ^/q  ^ according  to  (59), 

after  introducing  continuity,  reduce  to  (Y^/d)^.  The  ratios  according 

to  (57)  and  (58)  depend  on  d /d  and  U in  addition  to  Y./d.  Insofar 

x r C 1 

as  they  agree  more  closely  with  the  measured  values  of  q ^ /q  ^ for  Runs 
CS  and  CR  in  column  (8),  the  only  ones  for  which  there  were  equivalent 


U - Uc,  m/sec 

Figure  14.-  Variation  of  unit  total  sediment  discharge  with  effective 
velocity. 
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ice-covered  and  free-surface  flows,  the  estimates  based  on  (57)  and  (58) 
are  to  be  preferred.  From  these,  it  is  concluded  that  the  total  sedi- 
ment discharge  in  an  ice-covered  channel  may  be  reduced  to  less  than 
one  quarter  of  that  in  an  equivalent  free  surface  flow  if  the  cover  is 
smooth,  and  to  well  under  ten  percent  if  the  cover  is  rough. 

H.  Bed  Profiles.  Quantitative  measures  of  the  bed-form 
characteristics  were  obtained  from  a statistical  analysis  of  the  zero- 
crossing distances  and  amplitudes  of  the  digitized  longitudinal  bed 
profiles.  The  analysis  adopted  herein  is  essentially  as  described  by 
Annambhotla  (1969) . 

Longitudinal  profiles  taken  along  the  centerline  and  plotted 
by  the  Incremental  Plotter  at  the  University  Computer  Center  are  shown 
in  figures  16.  The  two  sets  of  profiles  for  the  A and  B series  of  runs 
were  taken  respectively  before  and  after  the  velocity  and  suspended 
sediment  profiles.  Insofar  as  the  profiles  appeared  to  indicate  a 
reasonably  stationary  stochastic  process,  they  were  not  passed  through 
a filter  to  attenuate  the  low  frequency  components  as  was  done  with 
Annambhotla' s Missouri  River  bed  profiles. 

A definition  sketch  for  the  zero-crossing  distances  and  ampli- 
tudes analysis  of  a bed  profile  is  shown  in  figure  17.  A zero-crossing 
is  defined  as  a point  where  the  bed  profile  crosses  the  line  defining 
the  mean  bed  level  in  a cartesian  coordinate  system  that  is  tilted  at 
the  flume  slope,  S.  In  figure  17  this  line  is  coincident  with  the  x-axis. 
The  zero-crossing  distance  is  then  defined  as  the  longitudinal  distance 
between  successive  zero-crossings.  A bed-form  length,  L,  is  the  sum  of 
two  successive  zero  crossings.  The  positive  and  nagative  amplitudes,  a+ 
and  a , are  defined  as  the  maximum  positive  and  negative  displacements 
of  the  bed  from  the  mean  level  between  two  successive  upward  (or  downward) 
zero  crossings.  The  height,  H,  of  a bed  form  is  the  sum  of  the  magni- 
tudes of  a pair  of  adjacent  positive  and  negative  amplitudes. 

Mean  values,  unbiased  standard  deviations,  and  unbiased 
skewness  coefficients  were  computed  for  H,  L,  and  a+,  and  |a  | from  each 
bed  profile.  For  a set  of  N discrete  values  Z1 , in  this  case 
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Longitudinal  bed  profiles. 
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Figure  16.  (continued). 
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representing  any  one  of  the  above  four  bed-form  characteristics,  the 
mean,  p,  is  computed  by 


W = N itfi 


the  unbiased  standard  deviation,  o,  by 


_ _ r 2 ,1/2 

0 = ta1 


and  the  unbiased  skewness  coefficient,  0,  by 


(N-l)(N-2)  2/3 

m2 

where  m2  and  m^  are  the  second  and  third  central  moments  defined 


respectively  by 


1 N 2 

"2  - if  jSi'y) 


1 N 3 

"3  * if  j5l(Z;fU) 


The  results  are  listed  in  table  8,  where  values  of  the  coefficient 
of  variation,  a/p,  are  listed  instead  of  the  standard  deviation.  For 
the  A and  B series  of  runs,  the  values  listed  are  the  averages  obtained 
from  6 profiles,  one  each  taken  along  the  centerline  and  the  sailing 
lines  one  quarter  of  the  flume  width  to  the  right  and  left,  before  and 
after  the  velocity  and  suspended  sediment  concentration  profiles  were 
obtained.  Comparison  of  the  results  for  the  individual  profiles  indi- 
cates no  consistent  significant  differences  between  those  for  the  cen- 
ter, left,  and  right  sailing  lines.  Therefore  profiles  were  taken  only 
along  the  centerline  in  the  C series.  Means,  standard  deviations,  and 
skewness  coefficients  of  H,  L,  a+,  and  |a  j for  each  profile  obtained 
in  the  A and  B series  are  tabulated  by  Song  (1978) . 
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Table  8.  STATISTICAL  MEASURES  OF  BED-FORM  CHARACTERISTICS 


Run 

Statistical 

Bed-form 

characteristics 

Measure 

Wave 

Wave 

Positive 

Negative 

height 

H 

length 

L 

amplitude 

a+ 

ampli tude 

Ui 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

AF 

Mean,  (m) 

0.0185 

0.216 

0.0086 

0.0100 

Coeff.  of 

Var. 

0.541 

0.627 

0.79 

0.82 

Skewness 

0.12 

0.72 

0.67 

0.71 

AS 

Mean,  (m) 

0.0184 

0 . 206 

0.0086 

0.0098 

Coeff.  of 

Var . 

0.565 

0.681 

0.79 

0.74 

Skewness 

0.32 

1.42 

0.76 

0.63 

AR 

Mean,  (m) 

0.0159 

0.200 

0.0073 

0.0086 

Coeff.  of 

Var. 

0.547 

0.680 

0.78 

0.  80 

Skewness 

0.21 

1.04 

0.62 

0.75 

BF 

Mean,  (m) 

0.0198 

0.312 

0.0088 

0.0110 

Coeff.  of 

Var. 

0.677 

0.839 

0.84 

0.92 

Skewness 

0.42 

1.25 

0.67 

0.85 

BS 

Mean,  (m) 

0.0201 

0.272 

0.0094 

0.0107 

Coeff.  of 

Var. 

0.562 

0.715 

0.74 

0.81 

Skewness 

0.19 

1.23 

0.54 

0.58 

BR 

Mean , (m) 

0.0176 

0.208 

0.0085 

0.0092 

Coeff.  of 

Var. 

0.574 

0.708 

0.76 

0.80 

Skewness 

0.21 

1.39 

0.55 

0.75 

CF 

Mean,  (m) 

0.0221 

0.388 

0.0097 

0.0124 

Coeff.  of 

Var. 

0.747 

0.750 

0.84 

1.10 

Skewness 

0.74 

0.73 

0.65 

1.22 

CS 

Mean,  (m) 

0.0208 

0.272 

0.0095 

0.0113 

Coeff.  of 

Var. 

0.587 

0.729 

0.74 

0.83 

Skewness 

0.22 

1.12 

0.36 

0.87 

CR 

Mean,  (m) 

0.0189 

0.218 

0.0091 

0.0098 

Coeff.  of 

Var. 

0.577 

0.726 

0.84 

0.84 

Skewness 

0.16 

1.26 

0.64 

0.81 

/ 
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The  following  trends  can  be  seen  from  the  results  in  table  8. 
Going  from  the  free-surface  (F)  to  the  smooth-cover  (S)  condition,  the 
mean  L,  the  coefficients  of  variation  for  H,  a+,  and  |a  |,  and  the  skew- 
ness coefficient  for  |a  | all  tend  to  decrease;  and  the  mean  H,  a+ , and 
|a  |,  and  the  coefficient  of  variationfor  L all  stay  about  the  same  (i.e. 
change  by  less  than  about  10  percent).  Going  from  the  smooth-cover  to 
the  rough-cover  condition,  the  mean  H,  |a  |,  and  L all  tend  to  decrease; 
and  the  coefficients  of  variation  for  H, |a  |,  and  L all  stay  about  the 
same.  The  overall  effect  of  the  simulated  ice  cover  appears  to  be 
toward  a reduction  and  uniformity  of  the  ripple  sizes.  Because  the 
flows  with  and  without  simulated  ice  covers  were  equivalent  only  in  the 
C series  of  runs,  the  tendencies  described  above  were  considered  to 
be  trends  only  if  they  occurred  in  the  C series  and  in  at  least  one  of 
either  the  A or  B series.  Whatever  effect  the  covers  may  have  had  on 
the  skewness  coefficients  appears  to  have  been  inconsistent  except  for 
the  effect  on  Ja  | as  noted  above.  However,  it  is  noteworthy  that 
all  skewness  coefficients  are  positive,  indicating  a greater  preponderance 
of  smaller  as  opposed  to  larger  ripples.  In  addition  there  is  a con- 
sistent tendency  for  the  skewness  coefficients  for  the  different  bed-form 
characteristics  to  increase  in  the  following  order:  H,  a+,  |a  |,  L. 

Percentage-distribution  histograms  for  H,  L,  a+,  and  |a  | for 
one  centerline  bed  profile  from  Runs  CF,  CS,  and  CR  are  shown  in  figure 
18.  No  attempt  was  made  to  fit  them  by  the  probability  density  functions 
corresponding  to  any  probability  law.  They  are  included  to  add  visual 
reinforcement  to  the  results  summarized  in  table  8.  Similar  histograms 
for  the  A and  B series  of  runs  are  shown  in  Song  (1978). 

The  relationship  between  the  bed-form  friction  factor  f"  and  the 
- -2 

relative  roughness  Y^L/H  based  on  the  mean  ripple  height  and  length  is 
shown  in  figure  19.  Experimental  results  from  the  present  investigation 
are  compared  there  with  data  for  free-surface  flows  obtained  in  the  same 
flume  and  with  the  same  sand  by  Annambhotla  (1969),  and  with  the  equation 


1 


3.3 


2.3 


(65) 


ibution  histograms  for  bed-form  characteristics.  Runs  CF , CS  and  CR 


(continued) . 


Positive  amplitude  distributions 


Negative  amplitude  distributions 


Annambhotla  (1969) 


ripple  dimensions. 
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Equation  65  is  a slightly  modified  form  of  Vanoni  and  Hwang's  (1967) 
logarithmic  bed-form  friction  factor  relation 


3.5  log  — - 2.3 
eH 


(66) 


where  e is  the  areal  concentration  of  ripples  on  the  bed.  The  changes 
are:  the  approximation  of  e by  H/L;  decreasing  the  value  of  the  numeri- 

cal coefficient  from  3.5  to  3.3;  and  the  substitution  of  Y,  for  r,  . Vanoni 

i b 

and  Hwang  found  that  the  equation  with  the  first  two  modifications 

provided  a rough  fit  to  the  data  from  several  investigations  wherein 

no  direct  measurements  to  determine  e had  been  obtained.  The  third 

modification  is  introduced  herein  as  another  test  of  the  hypothesis  that 

the  flow  in  the  lower  layer  of  an  ice-covered  channel  is  like  a free- 

surface  flow  with  depth  equal  to  Y^.  Except  for  the  two  points  of  Anna- 

bhotla's  that  are  far  to  the  left  of  the  others,  the  degree  of  scatter 

in  figure  19  does  not  much  exceed  that  in  Vanoni  and  Hwang's  original 

relationship.  It  is  therefore  concluded  that  Y^  is  about  as  good  a 

characteristic  flow  dimension  to  use  in  defining  the  relative  roughness 

of  the  bed  of  an  ice-covered  flow,  as  is  d or  r,  for  a free-surface  flow. 

b 


IV.  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

A.  Summary.  When  an  ice  cover  forms,  flow  processes  in 
alluvial  channels  assume  an  additional  degree  of  complexity.  The  im- 
mediate effect  of  the  ice  cover  is  to  approximately  double  the  wetted 
perimeter  of  the  channel,  which  radically  alters  the  distribut  Ions  of 
boundary  and  internal  shear  stresses,  and  of  velocity.  These  altera- 
tions in  turn  give  rise  to  several  additional  changes  which  usually 
include  an  increase  in  depth,  a reduction  in  velocity  and  sediment  trans- 
port, and  a change  in  the  bed  configuration  and  consequently  the  hydraulic 
roughness  of  the  bed.  The  extent  of  these  changes  evidently  depend  in 
large  measure  on  the  hydraulic  roughness  of  the  underside  of  the  ice 
cover  relative  to  the  roughness  of  the  bed.  However,  the  roughness  of 
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the  bed  depends  on  the  bed  configuration  which  in  turn  is  free  to  adjust 
to  the  new  shear  stress  distribution,  creating  an  internal  feedback  loop. 

The  main  objective  of  this  investigation  was  to  obtain  basic 
information  on  how  ice  covers  affect  the  sediment  transport  character- 
istics of  alluvial  streams  with  a view  toward  formulating  recommendations 
that  would  help  the  U.S.  Geological  Survey  to  improve  field  sampling 
and  data  collection  procedures  for  ice-covered  alluvial  channels.  To 
this  end,  a series  of  experiments  was  performed  in  a laboratory  flume 
wherein  detailed  comparisons  were  made  between  the  characteristics  of 
alluvial-channel  flow  processes  occurring  under  simulated  ice  covers,  and 
those  occurring  in  equivalent  free  surface  flows.  In  particular,  the 
experimental  investigation  addressed  the  effects  of  smooth  and  rough  ice 
covers  on  the  following  characteristics  of  the  flows:  (1)  depth  versus 
discharge  and  friction-factor  relationships;  (2)  velocity  distribution;  (3) 
suspended  sediment  concentration  profiles  and  discharge;  (4)  total  bed- 
material  discharge;  and  (5)  characteristics  of  bed  forms  and  their 
relationship  to  bed  friction  factors.  The  experiments  were  restricted 
to  flows  with  a 0.25  mm  sand  bed  in  the  ripple  regime. 

B.  Conclusions.  The  results  of  the  experiments  combined 
with  basic  principles  from  mechanics  and  hydraulics  gives  rise  to  the 
following  conclusions: 

1.  The  flow  in  an  ice-covered  channel  is  divided  by  a plane 
of  zero  shear  stress,  i.e.,  locus  of  maximum  velocity, 
into  a lower  and  an  upper  layer.  Since  the  shear  stress 
at  the  bed  of  an  ice-covered  channel  is  the  same  as  the 
shear  stress  at  the  bed  of  a free-surface  flow  with  the 
same  depth  and  average  velocity  as  the  lower  layer  of 
the  ice-covered  flow,  it  follows  that  features  of  the  flow 
which  depend  mainly  on  the  bed  shear  stress  and  the  mean 
velocity  should  also  be  essentially  the  same.  To  the 
extent  that  this  is  true,  it  should  be  possible  to  use 
relationships  developed  for  free-surface  flows  in  alluvial 
channels  for  predicting  certain  features  of  flows  in 
ice-covered  alluvial  channels. 
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2.  The  Karman-Prandtl  logarithmic  velocity  defect  law, 
commonly  used  for  free-surface  flows,  was  found  to  apply 
reasonably  well  to  both  the  lower  and  upper  layers  of 
flows  with  simulated  ice  covers.  The  vertical  distance 
from  the  lower  or  upper  boundary  to  the  locus  of  maximum 
velocity  is  simply  substituted  for  the  depth.  However, 
values  of  the  von  Harman  k may  differ  from  those  in  free- 
surface  flows.  Another  significant  feature  of  the  velocity 
distribution  observed  in  the  experiments  is  that  the  mean 
velocities  in  the  lower  and  upper  layers  are  nearly  equal 
to  each  other  and  to  the  overall  mean  velocity. 

3.  The  relative  thicknesses  of  the  lower  and  upper  layers, 
which  depend  mainly  on  the  roughnesses  of  the  lower  and 
upper  boundaries  relative  to  each  other,  were  found  tc 
depend  primarily  on  the  roughness  of  the  underside  of  the 
simulated  ice  cover  and  to  be  relatively  independent  of 
unit  discharge  over  the  range  of  experimental  conditions, 
implying  that  the  bed  roughness,  given  a particular  ice 
cover,  may  also  be  relatively  independent  of  unit  discharge. 

4.  The  total  depth  of  flow  in  a channel  with  a given  unit 
discharge  and  slope  was  found  to  be  significantly  increased 
by  the  addition  of  a simulated  ice  cover.  The  amount  of 
increase,  relative  to  the  depth  for  a free-surface  flow  con- 
dition, was  determined  to  be  from  20  to  30  percent  for  a 
smooth  cover  and  from  30  to  80  percent  for  a rough  cover. 

5.  The  Alam-Lovera-Kennedy  bed-form  friction  factor  relation- 
ships developed  for  free-surface,  alluvial-channel  flows, 
were  found  to  apply  equally  well  to  ice-covered,  alluvial- 
channel  flows  when  the  thickness  of  the  lower  layer  in 

an  ice-covered  flow  is  substituted  for  the  full  depth  of 
a corresponding  free-surface  flow  having  the  same  mean 
velocity  and  slope. 

6.  Solutions  of  the  Schmidt-0 ' Brien  equation,  with  turbulent 
mass  transfer  coefficients  for  suspended  sediment  based 

on  Bouissinesq  eddy  dif fusivities  computed  for  the  shear-stress 
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distribution  in  an  ice-covered  flow,  were  found  to  be 
qualitatively  compatible  with  the  measured  suspended 
sediment  concentration  profiles.  However,  good  quantita- 
tive agreement  between  computed  and  measured  profiles 
in  most  cases  could  only  be  obtained  by  utilizing 
unreasonably  large  values  of  the  ratio  of  sediment-to- 
momentum  diffusivity  (and/or  small  values  of  sediment 
settling  velocity)  in  the  computations.  Overall,  the 
results  of  the  suspended  sediment  concentration  measure- 
ments were  not  sufficiently  consistent  to  reach  any 
significant  conclusions  concerning  the  effect  of  an  ice 
cover  on  the  vertical  distribution  of  suspended  sediment. 

7.  The  total  bed-material  discharge  of  a flow  was  found  to 
be  greatly  reduced  by  the  addition  of  a simulated  ice 
cover.  Relative  to  total  sediment  discharges  in  free- 
surface  flows  with  the  same  unit  discharge,  slope, 

and  bed-material  characteristics , total  sediment  discharges 
were  determined  to  be  on  the  order  of  20  percent  for  flows 
with  a smooth  cover,  and  5 percent  for  flows  with  a rough 
cover. 

8.  The  overall  effect  on  bed  configuration  of  adding  a simu- 
lated ice  cover  was  a tendency  toward  a significant  re- 
duction in  ripple  sizes,  with  the  mean  length  decreasing 
somewhat  more  than  the  mean  height,  accompanied  by  a less 
pronounced  tendency  for  the  ripples  to  become  more  uniform 
in  size.  A rough  cover  tended  to  produce  greater  reduc- 
tion in  ripple  sizes  than  a smooth  cover,  but  the  uniformity 
of  ripple  sizes  was  about  the  same  for  smooth  and  rough 
covers . 

9.  The  thickness  of  the  lower  layer  was  found  to  be  about 

as  good  a characteristic  flow  dimension  to  use  in  defining 
the  relative  roughness  of  the  bed  of  an  ice-covered  flow, 
as  is  the  depth  or  hydraulic  radius  for  a free-surface  flow. 
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C.  Recommendations.  With  reference  to  possible  modifications 
in  the  field  sampling  and  data  collection  procedures  of  the  U.S. 

Geological  Survey  for  ice-covered  alluvial  channels,  the  results  of 
this  investigation  suggest  the  following: 

1.  Insofar  as  the  thickness  of  the  lower  layer  of  the  flow 
was  found  to  be  an  important  parameter  for  flows  with  an 
ice  cover,  it  is  recommended  that  velocity  measurements 
in  the  field  be  obtained  at  enough  additional  points  in 
selected-verticals  to  define  the  locus  of  maximum  velocity 
(i.e.  zero  shear  stress)  at  cross  sections  where  current- 
meter  discharge  measurements  are  obtained. 

2.  Any  information  on  bed  configuration,  even  if  only  visual 
(e.g.  whether  bed  appears  to  be  flat,  or  covered  by  dunes 
or  ripples)  would  be  useful  in  interpreting  data,  as  would 
be  descriptions  of  the  roughness  of  the  underside  of  the 
ice  cover. 

3.  Insofar  as  the  evidence  indicates  that  sediment  discharge 
is  substantially  reduced  by  ice  cover,  less  frequent  sedi- 
ment discharge  measurements  when  an  ice  cover  exists  would 
appear  to  be  justified  for  most  sampling  stations.  However, 
to  learn  more  about  sediment  transport  phenomena  in  ice- 
covered  streams,  fewer  sets  of  more  detailed  and  complete 
measurements  are  recommended. 
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ON  MEASUREMENT  OF  SUSPENDED  SEDIMENT  TRANSPORT 
IN  ICE-COVERED  ALLUVIAL  STREAMS, 

NOVEMBER  25,  1975 

The  Iowa  Institute  of  Hydraulic  Research  (IIHR)  was  commissioned  by  the 
USGS  to  review  and  evaluate  the  aspects  of  the  USGS  data  collection  and 
interpretation  program  that  relate  to  the  measurement  of  suspended-sedi- 
ment transport  in  ice-covered  alluvial  channels.  In  conjunction  with 
this  project  a one-day  workshop  was  held  at  The  University  of  Iowa,  Iowa 
City,  on  Tuesday,  November  25,  1975. 


Purpose.  The  main  purpose  of  the  workshop  is  to  exchange  information, 
identify  problems  and  possible  solutions  thereto,  and  formulate  goals 
on  matters  relating  to  sampling  techniques  and  equipment,  data  analysis 
and  interpretation,  and  design  of  data  collection  programs.  The  work- 
shop will  serve  the  dual  purpose  of  providing  IIHR  project  personnel 
with  sources  of  information  concerning  USGS  field  experience  in  ice-covered 
alluvial  streams,  and  providing  USGS  personnel  with  information  and  in- 
sights resulting  from  IIHR  investigations  of  ice  regimes  and  the  mech- 
anics of  flow  in  ice-covered  channels. 

List  of  Participants 

John  Ficke,  USGS,  Reston,  Virginia 
Steve  Hindall,  USGS,  Madison,  Wisconsin 
Carl  Nordin,  USGS,  Denver,  Colorado 
John  Skinner,  USGS,  Minneapolis,  Minnesota 
Jack  Kennedy,  IIHR,  University  of  Iowa 
Bill  Sayre,  IIHR,  University  of  Iowa 
Ed  Slattery,  IIHR,  University  of  Iowa 
George  Ashton,  CRREL,  Hanover,  New  Hampshire 

Agenda.  The  meeting  will  be  loosely  structured  to  allow  for  free 
interchange  and  spontaneity.  A tentative  list  of  the  topics  to  be  cov- 
ered follows: 

1.  Sediment  transport  in  ice-covered  alluvial  streams  within  context  of 
USGS  water  measurement  and  sampling  program. 
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2.  Measurement  of  suspended  sediment  transport  in  ice-covered  channels. 

a.  equipment  and  procedural  problems  due  to  ice  cover  and/or  freezing 
conditions. 

b.  presentation  of  typical  sets  of  data  that  illustrate  the  effects 
of  ice  covers  on  sediment  concentration  and  velocity  distribution 
data. 

c.  design  of  sampling  program,  location  of  sampling  verticals  and/or 
points,  frequency  of  sampling. 

d.  analysis  and  Interpretation  of  data. 

3.  Measurement  of  supporting  data  such  as  discharge,  slope,  bed  con- 
figuration, bed  load  transport,  ice  ripples,  ice  texture.  How  much 
of  this  kind  of  supporting  data  is  needed  to  satisfy  USGS  objectives? 

4.  Mechanics  of  flow  in  ice-covered  channels. 

a.  ice  regimes — formation,  evolution,  and  breakup  of  ice  covers. 

b.  resistance  to  flow,  distribution  of  velocity  and  suspended  sedi- 
ment under  different  ice  regimes. 

c.  speculations  on  Interactions  involving  ice  cover,  bed  configuration, 
and  bed  and  suspended  load  transport. 

5.  Tour  of  IIHR  laboratories  with  emphasis  on  ice  and  sediment  trans- 
port facilities. 

Assignments.  Everyone  will  no  doubt  have  something  to  contribute  on  each 
agenda  item.  However,  the  following  assignments  are  being  made  to  ensure 
that  all  bases  are  covered  and  to  give  each  participant  a clearer  concept 
of  his  role  and  the  particular  contributions  he  is  expected  to  make. 

Flcke: 

1.  Introductory  comments  relating  to  item  1 on  agenda. 

2.  Contribute  perspective  to  items  2c,  2d  and  3 from  viewpoint  of  total 
program  objectives. 

Hindall : 

1.  Contribute  benefit  of  field  experience  in  measuring  sediment  transport 
and  related  phenomena  in  ice-covered  channels,  and  experience  in 
analyzing  and  interpreting  the  data.  Should  be  prepared  to  make  a 20-  to 
30-mlnute  oral  presentation  with  slides,  photographs,  graphs,  tables, 
etc.  on  items  2a  and  2b  of  agenda,  to  get  the  discussion  started. 
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2.  Contribute  perspective  to  items  2c,  2d  and  3 from  viewpoint  of 
field  investigator. 

Ashton: 

1.  Oral  presentation  with  slides  (30  minutes)  covering  spectrum  of 
ice-related  phenomena  in  rivers. 

2.  Contribute  broad  perspective  on  river-ice  phenomena. 

Nordin: 

1.  Resource  person  on  relationships  to  other  USGS  alluvial  channel 
investigations . 

2.  Contribute  perspective  to  all  items. 

Skinner: 

1.  Resource  person  on  sampling  equipment  and  instrumentation. 

Kennedy  and  Sayre: 

1.  Oral  presentations  to  initiate  discussion  on  item  4 of  the  agenda. 

2.  Contribute  perspective  on  items  2c,  2d  and  3 from  the  academic  viewpoint. 

Slattery: 

1.  Audiovisual  equipment.  Transparency,  opaque  and  slide  projectors. 

Tape  recorder. 

2.  Contribute  perspective  from  the  viewpoint  of  the  graduate  research 
assistant  who's  going  to  be  doing  much  of  the  work  on  this  project. 

Travel  and  Accommodations.  The  most  convenient  way  to  reach  Iowa  City 
is  by  air  to  Cedar  Rapids.  Iowa  City  is  about  20  miles  from  the  Cedar 
Rapids  Airport.  Reservations  for  the  entire  group  have  been  made  for 
the  nights  of  November  24  and  25  and  the  Iowa  House  in  the  Iowa  Memorial 
Union,  located  on  the  University  campus.  Direct  transportation  between  the 
Airport  and  the  Iowa  House  is  provided  by  Iowa  Limousine  Service,  Inc. 
for  all  incoming  and  outgoing  flights,  so  it  is  not  necessary  to  rent  cars 
at  the  Airport. 

Social  Activities.  All  of  the  workshop  participants  are  invited  to  dinner 
on  the  evening  of  Monday,  November  24,  at  one  of  the  well  known  Amana  Colony 
restaurants,  as  guests  of  the  Iowa  Institute  of  Hydraulic  Research.  We 
will  meet  in  the  lobby  of  the  Iowa  House  from  where  transportation  will 
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be  provided,  at  7:00  p.m.  Please  note  when  making  travel  reservations. 

lime  and  Place  of  Workshop.  8:00  a.m. , Tuesday,  November  25,  in  Michigan 
State  Room,  third  floor  of  Iowa  Memorial  Union. 


Information.  For  additional  information,  if  needed, 
Sayre  (319-353-3082)  or  Jack  Kennedy  (319-353-4679). 


call  either  Bill 
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FIELD  TRIP  TO  UPPER  IOWA  RIVER  NEAR  DORCHESTER,  IOWA 
JANUARY  20,  1976 

Personnel:  Bill  Matthes  and  Dick  Downs,  USGS, 

Iowa  District  Office,  performed  measurements. 

Bill  Sayre  and  Ed  Slattery,  IIHR, 

accompanied  as  observers  and  assisted  In  measurements. 

I.  Purpose 

A.  Observe  discharge  and  sediment  measurement  techniques  for  Ice- 
covered  streams. 

B.  Become  familiar  with  measurement  and  sampling  equipment. 

II.  Description  of  Gaging  Station 

A.  Located  at  Route  76  bridge 

B.  Gage  house  contains  stage  recorder  consisting  of  bubble  gage,  digital 
punched-tape  recorder,  and  strip-chart  recorder  for  backup. 

C.  Manual  plumb  bob  for  measuring  stage,  attached  to  wire  rope  and 
reel  assembly  mounted  on  bridge  railing. 

III.  Selection  and  Preparation  of  Cross  Section 

A.  Cross  Section  located  in  straight  reach  approximately  200  ft  up- 
stream bridge  where  streamlines  parallel  to  each  other  and  to 
banks,  and  composition  of  bed  considered  to  be  typical. 

B.  Probe  for  weak  spots  across  section  with  ice  chisel.  Stretch 
tagline  across  river  and  fix,  perpendicular  to  shore. 

C.  Locate  positions  of  sampling  verticals.  Drill  holes  with  power 
auger  along  tagline  at  about  2-ft  intervals  above  channel  where 
water  flowing,  and  at  about  4-  to  7-ft  intervals  in  shallow  areas 
where  ice  reaches  to  bottom  of  channel. 

IV.  Depth  Soundings 

A.  Measure  total  depth  from  top  of  ice  to  bottom  of  channel  at 
each  hole  with  ice  rod,  similar  to  surveyor's  rod. 

B.  Measure  thickness  of  ice,  using  steel  angle  attached  to  ice  rod  to 
find  underside  of  ice  cover. 

C.  Subtract  ice  thickness  from  total  depth  to  get  water  depth. 

V.  Velocity  Measurements 

A.  Velocities  measured  with  vane-type  current  meter  attached  to 
bottom  of  rod  supported  from  above  by  tripod  assembly  resting  on 
ice. 

B.  Velocities  measured  at  0.6  depth  from  underside  of  ice  where  water 
depth  equal  to  or  less  than  2.5  ft,  and  at  0.2  and  0.8  depth  where 
depth  greater  than  2.5  ft. 

C.  Beeps  counted  for  at  least  40  seconds. 
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0.  Velocities  measured  at  0.6  depth  multiplied  by  empirical  correc- 
tion factor  of  0.92  to  get  mean  velocity,  to  account  for  effect 
of  ice  cover  on  velocity  distribution. 

VI.  Suspended  Sediment  Sampling 

A.  Samples  obtained  at  about  8 verticals  using  a US  DH-48  hand 
sampler  mounted  to  a rod. 

B.  Depth-integrated  samples  obtained.  Except  at  deepest  sampling 
locations,  two  locations  sampled  with  one  bottle. 

C.  Very  little  suspended  sediment  in  samples. 

VII.  Bed-Material  Sampling 

A.  Piston-type  core  sampler  used  to  obtain  bed-material  samples 
at  same  verticals  where  suspended  sediment  samples  obtained. 

B.  Samples  placed  in  labeled  plastic  frozen-food  containers. 

C.  Samples  appeared  to  be  well-graded  with  significant  quantity  of 
coarse  sand. 

VIII.  Bed-Load  Sampling. 

A.  Helley-Smith  bed-load  sampler  used  for  measuring  bed  load  at  same 
sampling  locations  where  suspended  sediment  and  bed-material  samples 
obtained. 

B.  Sampling  holes  in  ice  had  to  be  enlarged  with  chain  saw  before 
sampler  could  be  lowered  into  water. 

C.  Sampler  lowered  to  bed  in  vertical  position,  left  on  bed  in 
horizontal  position  pointing  upstream  for  30  seconds,  and  raised 
to  surface  in  vertical  position. 

D.  Sample  transferred  from  collector  to  labeled  container. 

IX.  Ice  and  Bed  Features 

A.  Blocks  of  ice,  about  4 ft,  x 1 ft,  were  cut  with  chain  saw  and 
turned  over. 

B.  Observed  ripple-like  features  with  well-defined  ridge-like  crests 
on  undersurface  of  ice. 

C.  Water  clear  enough  to  see  bed  through  holes.  Bed  mostly  flat. 

Ripple  crests  visible  where  bed  composed  of  fine  sand. 

X.  Comments  and  Criticisms. 

A.  V-D.  Correction  factor  of  0.92  to  account  for  ice-cover  effect 
on  velocity  distribution  seems  very  arbitrary.  Shape  of  velocity 
profile  depends  on  roughness  of  ice  cover  relative  to  roughness 
of  bed,  which  can  vary  considerably. 

B.  VII-A.  Much  of  bed-material  sample  lost  while  raising  piston 
sampler.  More  difficult  to  recover  sample  without  losing  some  than 
under  free-surface  conditions  because  small  hole  in  ice  restricts 
freedom  of  motion  in  raising  sampler. 
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